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Abstract: 

The exact solution of the noncompact SL{2, C) Heisenberg spin magnet reveals a hidden sym- 
metry of the energy spectrum. To understand its origin, we solve the spectral problem for 
the model within quasiclassical approach. In this approach, the integrals of motion satisfy the 
Bohr-Sommerfeld quantization conditions imposed on the orbits of classical motion. In the rep- 
resentation of the separated coordinates, the latter wrap around a Riemann surface defined by 
the spectral curve of the model. A novel feature of the obtained quantization conditions is that 
they involve both the a— and /3— periods of the action differential on the Riemann surface, thus 
allowing us to find their solutions by exploring the full modular group of the spectral curve. 
We demonstrate that the quasiclassical energy spectrum is in a good agreement with the exact 
results. 
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1. Introduction 



Exact solution of the spectral problem for quantum-mechanical multi-particle systems is the 
central problem in the theory of Integrable Lattice Models [|^ . One of the best known examples 
of such systems is spin— 1/2 Heisenberg spin magnet. The model can be solved exactly by the 
Algebraic Bethe Ansatz (ABA) and it has numerous applications |0, ^ The Heisenberg 
magnet model can be generalized from the compact SU (2) spins to noncompact spins, "living" 
in an infinite-dimensional, unitary representations of the SL[2, C) group. The corresponding 
integrable model describes a nearest-neighbour interaction between the SL{2, C) spins and is 
called the noncompact Heisenberg magnet [§. 

The noncompact Heisenberg spin magnets have important implications in high-energy QCD [^, 
0]. It is well-known that hadronic scattering amplitudes grow as a power of energy in agree- 
ment with the Regge model. In perturbative QCD framework, this behaviour can be attributed 
to a contribution of colour-singlet gluonic compound states. These states satisfy the Bartels- 
Kwiecinski-Praszalowicz equation which coincides, in the multi-colour limit, with the Schrodinger 
equation for the noncompact SL{2, C) spin magnet. The effective QCD interaction between N 
reggeized gluons (A^ = 2, 3, . . .) occurs on two-dimensional plane of transverse degrees of freedom 
(the impact parameter space). It is described by the Hamiltonian Hn, which defines a quantum- 
mechanical system of N particles with the coordinates Zk = {xk, Uk) and the momenta pk = —idk, 
such that [zk,Pn] = i^kn^"^ {k,n = 1,. . . ,N and a,P = 1,2) with the Planck constant h = 1. 
To map this system into a Heisenberg magnet, one introduces holomorphic, z = x + iy, and 
antiholomorphic, z = x — iy, complex coordinates on the plane and defines the spin operators as 

Sk = izkPk + s, S^=-ipk, = izlpk + 2szk . (1.1) 

Here pk = —id/dzk is a (complex) momentum along the 2;— direction and the parameter s = 
(1 + ns)/2 + iug (with Ug integer and real) is a single-particle SL{2, C) spin. One also defines 
the antiholomorphic spin operators 5°, 5*^ and 5^ acting along the z— direction. They are given 
by similar expressions with Zk and s replaced by Zk = zl and s = 1 — s* = (1 — ns)/2 + iug, 

respectively. The operators Sk and Sn act along different directions on the z— plane and therefore 
commute. They satisfy the standard sl{2) commutation relations 

[^k^ ] = '^Sl^kn , [5*°, S^] = ±S^6kn , (1-2) 

with the quadratic Casimir operator S*^ = (5°)^ + (S+S*" + S^S^)/2 = s{s — l). Similar relations 

hold in the antiholomorphic sector. Defined in this way, the spin operators Sk and Sk are the 
generators of the unitary principal series representation of the SL{2, C) group labelled by the 
pair of spins {s, s), or equivalently integer rig and real Ug 0. The interaction between N reggeized 
gluons in multi-colour QCD is translated into the nearest-neighbour interaction between the spins 
Sk and Sk for the SL{2, C) Heisenberg magnet at s = and s = 1 ^ 

N 

'Hn = Hk,k+i 1 Hk,k+i = H{Sk ■ Sk+i) + H{Sk ■ Sk+i) , (1-3) 

k=l 

where Hn,n+i = Hn,i and the two-particle Hamiltonian is expressed in terms of the Euler 
-0— function, H{x) = ip{j{x)) + ip^l — j{x)) — 2'i/'(l) with j(j — 1) = 2a; + 2s{s — 1). As follows 
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from (|1.3|) , TIn can be split into a sum of two mutually commuting Hamiltonians acting along 
the z— and z— directions. 

The noncompact SL{2, C) Heisenberg magnet ( |1.3D is a completely integrable model. It 
possesses a large enough set of mutually commuting conserved charges g„ and g„ (n = 2, A^) 
such that g„ = and [7iAr,g„] = [HNyQn] = 0. The charges g„ are polynomials of degree n in 
the holomorphic spin operators. They have a particular simple form at s = P, 0] 

In = ^jlj2^j2j3---^jnjlPjlPj2---Pjn i^'^) 

l<jl<j2<...<jn<N 

with = -2j — and defined in ([1.11). The "lowest" charge q2 is related to the total spin of 
the system h. For the principal series of the SL{2, C) it takes the following values ^ 

q2 = -h{h-l)+Ns{s-l), h = ^-^ + iUh, (1.5) 

with Uh integer and Uh real. The eigenvalues of the integrals of motion, q2,...,qN, form the 
total set of quantum numbers parameterizing the eigenstates of the model (|1.3|) . As was already 
mentioned, at s = and s = 1 the latter define the multi-gluonic compound states in multi-colour 
QCD. 

In spite of the fact that the noncompact SL{2, C) Heisenberg magnet represents a generaliza- 
tion of the compact SU{2) spin chain, a very little has been known about its energy spectrum till 
recently. One of the reasons is that the exact solution of the eigenproblem for the Hamiltonian 
( |1.3|) represents a difficulty of principle. In distinction with the compact magnets, the quantum 
space of the SL{2, C) magnet does not possess the highest weight and, as a consequence, the 
conventional methods like the ABA method 0, ^ are not applicable. 

The eigenproblem for the noncompact SL{2, C) Heisenberg magnet has been solved exactly in 
Refs. 1^, 1^, |l^] using the method the Baxter Q— operator ||l|]. This method allowed us to establish 
the quantization conditions for the integrals of motion of the model, gs, ^at, obtain an explicit 
form of the dependence of the energy on the integrals of motion, Ejy = E^ivh^ n^, gs, ^at), and 
construct the corresponding eigenfunctions in the representation of the Separated Variables ||11|| . 
Solving the quantization conditions, we calculated the spectrum of the noncompact SL{2, C) 
magnet of spin s = for the number of particles 2 < A^ < 8. Its close examination revealed the 



following properties of the model |TD 



Quantized values of the charges qk (with k = 3,...,N) depend on the "hidden" set of 
integers £ = (£i, £2, ^2{Ar-2)) 

qk = qki^h;nh,e) , (1.6) 
where integer rih and real Uh define the total SL{2, C) spin of the state, Eq. ( |1.5| ). 



As a function of continuous Uh, the charges form the family of trajectories in the moduli 
space q = {q2,<l3, ■■■,<1n) labelled by integers and £ . Each trajectory in the g— space 
induces the corresponding trajectory for the energy E]\r (see Figure below) 

En = EN{vh\nh,t) . (1.7) 

For fixed total SL{2, C) spin of the model, Eq. ( |1.5|) , the eigenvalues of the "highest" charge 
g^^ define an infinite set of distinct points on the complex g^^— plane. At A^ = 3 and 
A^ = 4 they are located close to the vertices of a lattice built from equilateral triangles and 
squares, respectively (see Figures H and || below). 
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Their origin remains obscure mainly due to a rather comphcated form of the exact quantization 
conditions. The main goal of this paper is to present a physical interpretation of these properties. 
Our analysis is based on a generalization of the well-known quasiclassical methods to noncompact 
Heisenberg magnets. One might expect a priori that these methods could be applicable only for 
high excited states. Nevertheless, as we will demonstrate below, the quasiclassical formulae work 
with a good accuracy throughout the whole spectrum of the noncompact SL{2, C) Heisenberg 
magnet. 

To formulate the quasiclassical solution of the eigenproblem for Hamiltonian (|1.3| ), one has 
to consider a classical analog of the noncompact Heisenberg spin magnet [|I^, |TE[. From point 
of view of classical dynamics, the model describes a chain of interacting particles on the two- 
dimensional i*— plane. We use (anti)holomorphic variables on the phase space and define the 
coordinates and the momenta of particles as Zk = {zk, Zk) and pk = {pk,Pk), respectively. By the 
definition, Zk and pk take complex values such that Zk = zl and pk = p^.Q The only non-trivial 
Poisson bracket is given by {zk,Pn} = {zk,Pn} = ^kn- The classical model inherits a complete 
integrability of the quantum noncompact spin magnet. Its Hamiltonian and the integrals of 
motion are obtained from ( |1 . 1| ) and ( |1.4[ ) by replacing the momentum operators by the 

corresponding classical functions leading to {g^, TIn} = {%, <in} = 0. Since the Hamiltonian (|1.3|) 
is given by the sum of holomorphic and antiholomorphic functions, from point of view of classical 
dynamics the model describes two copies of one- dimensional systems "living" on the complex z— 
and ^— lines. The solutions to the classical equations of motion have a rich structure and turn 



out to be intrinsically related to the finite-gap solutions to the soliton equations [H, 15|. Namely, 
the classical trajectories have the form of soliton waves propagating in the chain of particles. 
Their explicit form in terms of the Riemann 6'— functions was established in [|1^] by the methods 
of the finite-gap theory |]14|, |T^. The charges q define the moduli of the soliton solutions and 



take arbitrary complex values in the classical model. Going over to the quantum model, one 
finds that q are quantized. In the quasiclassical approach presented in this paper, their values 
satisfy the Bohr-Sommerfeld quantization conditions imposed on the orbits of classical motion 
of particles. 

In a standard manner, the WKB ansatz for the eigenfunction of the model ( |1.3|) involves 
the "action" function, \1/wkb(^i, • • • , ^tv) ~ exp{iSo/h). Due to complete integrability of the 
classical system, it can be defined as a simultaneous solution to the system of the Hamilton- 
Jacobi equations 

^ dSo ^ f dSo 



k=l 



dzk ^' ^4^'^)=^"' in = 2,...,N), {U 



where z = {zi, zn) denotes the set of holomorphic coordinates, q„(2,p) stands for the symbol 
of the operator ( p..4|) and P is a holomorphic component of the total momentum of particles. 
The z— dependence of 5*0 is constrained by similar relations in the antiholomorphic sector. To 
find a general solution to Eq. (|1.8| ), one performs a canonical transformation to the classical 
separated coordinates |TI|, [I^ 

{zi,Z2, ...,Zn) {zo,Xi,X2, ...,xn-i) , (1-9) 



^Of course, one can work instead with real, Cartesian coordinates, but our choice is advantageous as it is 



dictated by the chiral structure of the Hamiltonian (1.3). 
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with zo the center-of-mass coordinate of the system and Xn = {xn,Xn = a;*) new collective 
(separated) coordinates. As explained in Section 2.1, the classical dynamics in the separated 
variables is determined by the spectral curve ( "equal energy" condition) 

r^v: = t%{x) - 4x^^ , t^ix) = 2x^ + q2X^-^ + ... + qn-iX + , (I-IO) 

with y{x) = 2x^ sinhpx and Px being the momentum in the separated coordinates. Here t^ix) 
is a polynomial of degree with the coefficients defined by the holomorphic integrals of motion 
Qn. The spectral curve establishes the relation between the holomorphic components of the 
separated coordinates, x and p^, for a given set of the "energies" q2, . . ., Qn- As we will see 
below, its properties play the central role in our analysis. 

In the separated coordinates, the solution to the Hamilton- Jacobi equations ( |1.8| ) takes the 
form So{zo, xi, X2, ...,xn-~i) = {P ■ zq) + Y.k=\ 'S'o(^fc) with lH 



5*0 (x) = / dxpx + / dxpx = 2'Re / dxp^. (l-H) 



Here complex momentum px = was defined in ( |1.10|) and xq is arbitrary. The WKB expression 



for the wave function in the separated coordinates, ~ exp {iSo{zo,xi,X2, ...,xn-i)/K), factorizes 
into a product of single-particle wave functions, QwrbI^^) ~ exp (2S'o(xfc)//i). According to 
( |1.10|) , the momentum, p^, and, as a consequence, the action function So{x) are multi- valued 
functions of x. Denoting the different branches of the action function as 5*0, a (x), one writes the 



WKB expression for the wave function of the quantum spin magnet as a sum over branches ||16|, p7? 



Qwkb(x) = ^ Aa{x) exp ( ^So^a{x) j . (1.12) 

The function Ak{x) takes into account subleading WKB corrections and is uniquely fixed by 
So,a{x). In general, the expression in the r.h.s. of ( |1.12| ) is not a single- valued function of x. 
For Qwkb{x) to be well-defined, the charges q have to satisfy the Bohr- Sommerf eld quantization 
conditions. One of the main results in this paper is that these conditions can be expressed in 
terms of the periods of the "action" differential over the canonical set of the a— and /3— cycles 
on the Riemann surface corresponding to the complex curve ( |1.10 ) 



Re® dxpx = 7Thi2k-i , Red) dxpx = Tihl2k , {k = 1, .., N — 2) , (1-13) 

with i = {ii, . . . , i2N~i) being the set of integers. 

The relations ( p..l3|) define the system of 2(A^ — 2) real equations on the (A^ — 2) complex 
charges q3,...,qN (we recall that the eigenvalues of the "lowest" charge q2 are given by ( |1.5| )). 
Their solutions lead to the quasiclassical expressions for the eigenvalues of the integrals of motion 
of the noncompact spin magnet. As we will demonstrate in Section 5, these expressions have the 
form ( |1.6| ) and are in a good agreement with the exact results of Ref. [^, |10|- A novel feature of the 
quantization conditions ( |1.13|) is that they involve both the a— and /9— periods on the Riemann 
surface. This should be compared with the situation in one- dimensional lattice integrable models. 



like the Toda chain model [16, [T^ and the 5'L(2,M) Heisenberg spin magnet [19, |2y, 21[. There, 
the WKB quantization conditions involve only the a— cycles, since the /3— cycles correspond to 
classically forbidden zones. In the SL{2, C) spin magnet, the classical trajectories wrap over an 
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arbitrary closed contour on the spectral curve ( |1.10|) leading to ( |1.13|) . This fact allows one to 
explore the full modular group ||2^ of the complex curve (|1.10| ). 

The paper is organized as follows. In Section 2 we remind the definition of the noncompact 
Heisenberg spin magnet both in the classical and quantum cases. Going over to the represen- 
tation of the Separated Variables we construct the wave function of the model in terms of the 
solutions to the Baxter equation. Applying the WKB methods, we solve the Baxter equation 
in Section 3 and show that the quasiclassical expression for the wave function is uniquely de- 
fined by the complex curve introduced in Requiring this function to be single- valued, 
we obtain the quantization conditions ( |1.13D for the integrals of motion In Section 4 we 
obtain quasiclassical expressions for the energy and the quasimomentum. Both observables are 
expressed in terms of the Q— blocks, which satisfy the holomorphic Baxter equation and have 
prescribed analytical properties and asymptotic behaviour at infinity. In Section 5 we analyze the 
quantization conditions (|1.13| ) and compare their solutions with the exact results for the energy 
spectrum. Section 6 contains concluding remarks. Some technical details of the calculations are 
summarized in two Appendices. 



2. Noncompact Heisenberg spin magnet 



Let us summarize, following |jT2|, |T3[, the main features of the SL(2, C) Heisenberg spin magnet 



in the classical and quantum mechanics. 



2.1. Classical model 



In the classical case, the model describes the chain of N interacting particles on the two- 
dimensional plane with the Hamiltonian (|1.3|) . The classical motion along the complex 2;— direction 
is described by the Hamilton equations 

dtZk = {zk, Hn} = , dtPk = {Pk, T-Cn} = -^J^ . (2.1) 

dpk dzk 

This system is completely integrable and the integrals of the motion are given by ( |1.4| ) . Following 
the Quantum Inverse Scattering Method , one can describe the classical Heisenberg spin magnet 
by the Lax matrix 

Lk{u) =u-t + iSl-a' + iSt ■ a + zS, ■ a+ = ( " + f ° ^ ^_^^^o ) , (2-2) 

with a"" being the Pauli matrices. The dynamical variables 5°, 5*^ and depend on the 
(holomorphic) coordinates and momenta of particles, Zk and pk, respectively, and they are given 
by the same expressions as in ( p..l| ). The Hamilton equations ( p.l|) are equivalent to the matrix 
Lax pair relation 

dtLk{u) = {Lk{u),HN} = Ak+i{u)Lk{u) - Lk{u)Ak{u) , (2.3) 

with Ak{u) being a 2 x 2 matrix depending on the coordinates and momenta of particles. 

The exact integration of the classical equations of motion is based on the Baker-Akhiezer 
function \E'fc(M;t) [|l^]. By the definition, it satisfies the system of matrix relations 

Lkiu)^kiu; t) = ^k+iiu] t) , dt^kiu] t) = Ak{u)'^k{u] t) . (2.4) 
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Introducing the monodromy matrix as a consecutive product of the Lax matrices, one finds that 
it produces the shift of the Baker- Akhiezer function along the chain 



Tn{u) = Ln{u)...Li{u), TN{u)^i{u;t) = ^N+i{u;t). (2.5) 

For periodic boundary conditions, Zk+N = Zk and pk+N = Pk, the Baker- Akhiezer function satisfies 
the Bloch-Floquet relation 

^N+kiu;t)=w{u)^kiu;t), (2.6) 

According to (|2.5| ), the Bloch-Floquet factor w{u) is an eigenvalue of the monodromy matrix. 
Therefore, it does not depend on the time and satisfies the characteristic equation 

det {Tn{u) -w) = w^ -wtN{u) + = . (2.7) 

Here t^iu) = tr T/v(u) is a polynomial in u of degree with the coefficients given by the integrals 
of motion, Eq. ( l.lOp . Introducing the complex function y{u) = w — u^^ /w, one obtains from 



that y{u) defines the algebraic complex curve (|1.10|) . 



The Baker- Akhiezer function "^kiu; t) is a double- valued function on the complex m— plane |T5 
Its explicit expression in terms of the theta-functions defined on the curve ( p..lO| ) can be found 
in |T^. We do not present it here since we will not use the Baker- Akhiezer function in the rest 
of the paper. The function \E'fc(M; t) has — 1 simple poles at u = Xk {k = 1, . . . , N — 1) and the 
same number of zeros. Remarkable property of its poles is that the variables {xk,p{xk)) (with 
p{u) = \n{w{u)/u'^) and A; = 1, . . . , A^ — 1) form the set of holomorphic separated variables for 
the classical model, Eq. ( |1.9|) . Notice that Xk and p{xk) take arbitrary complex values. This 



allows one to integrate the equations of motion exactly and reconstruct the classical trajectories 
of particles on the Riemann surface defined by the curve ( 1.10| ). The same classical motion 



describes a soliton wave propagating in the chain of A^ particles on the two-dimensional i*— plane. 
2.2. Quantum model 

In the quantum case, the eigenfunction of the SL{2, C) Heisenberg spin magnet, ^'(ii, . . . , zn), 
is defined as a simultaneous eigenstate of the integrals of motion q2,...,qN, Eq. ( p..4|) , and 
their antiholomorphic counterparts. Together with the total momentum of the system, P, their 
eigenvalues q = (g2, • • • , Qn) define the total set of the quantum numbers of the model. Due 
to chiral structure of the Hamiltonian and the integrals of motion, the eigenfunction can be 



decomposed as ||23| 



M/(fi, . . . , zn) = J2 Cabiq, g)^^'^^ {zi, . . . , Ziv)^f {zi, ...,Zm), (2.8) 

a, b 

where '^q'^^ and \E'^^^ diagonalize the integrals of motion in the holomorphic and antiholomorphic 
sectors, respectively, and Cab are mixing coefficients. The function "^{zi, . . . , z^) is a single- valued 
function on the two-dimensional plane. It belongs to the principal series of the SL{2, C) group 
labelled by the spins [h, h) defined in ( |1.5| ) (with h = 1 — h*) and is normalizable with respect 

to the SL{2, C) scalar product. In contrast with \E'(ii, . . . , zn), the chiral solutions "^^q^ and "^^^ 
acquire nontrivial monodromy when Zk encircles other particles on the plane. The quantization 
conditions for the integrals of motion q follow from the requirement that the monodromy should 
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cancel in the r.h.s. of ( p.8|) . The same condition fixes (up to an overall normalization) the mixing 
coefficients Cab- 
To formulate the quantization conditions it is convenient to switch from the coordinate 
^—representation to the representation of the Separated Variables (SoV) |[ll|]. In this repre- 
sentation, the wave function takes a factorized form 

^{Zi, Z2, zn) ^ $(io, Xi, X2, xn-i) = Q{xi) Q{X2) ■ ■ ■ Q{xn) , (2.9) 

where P is the total momenta of particles, zq is the center-of-mass coordinate of the system 
and xi,...,XAr_i are new collective (separated) coordinates. The exphcit form of the unitary 
transformation to the SoV representation can be found in (see also |^ for similar expressions 
at = 2 and = 3). The eigenfunction in the SoV representation has the following properties. 
Introducing holomorphic and antiholomorphic components x = [x, x), one finds that the possible 
values of the separated coordinates can be parameterized by integer n and real z/ as 

x = u- — , x = u+ — , (2.10) 

so that X = X* and i{x — x) = n. Here, as before, one has h = 1. To restore the ^—dependence 
one has to substitute n — > hn. Notice that, in contrast with the classical case, the separated 
variables have a discrete imaginary part for finite h. 

A single-particle Q— function entering ( p.9| ) satisfies the holomorphic Baxter equation 

{x + is)^Q{x + i, x) + {x — is)'^Q{x — i,x) = tN{x) Q{x, x) , (2.11) 

with tN{x) defined in (|1.10| ). Similar equation holds for Q{x,x) in the antiholomorphic sector 
with s and g„ replaced by s = 1 — s* and g„ = q*, respectively. The solution to the Baxter 
equations, (^(z/ — in/2, v + in/2), is a well-defined, regular function of integer n and real v. At 
large u and fixed n it has the following asymptotic behaviour 

Q{u - in/2, V + m/2) ~ e'® ^h^h-n^s^s) ^ ^-io ^i-h+i-h-Nis+s) ^ ^2.12) 

where h and h = 1 — h* define the total SL{2,C) spin of the model, Eq. (|1.5|), and is some 
phase. 

The exact solution to Eqs. (|2.11| ) and ( p.l2| ) was constructed in Refs. |T0[. It allowed 
us to establish the quantization conditions for the integrals of motion and calculate the energy 
spectrum of the model. In this paper we shall present another, quasiclassical approach to solving 
the Baxter equations. Although it does not provide the exact solution for the Q— function, it 
allows us to elucidate a hidden symmetry of the energy spectrum. 



3. Quasiclassical wave function 

The quasiclassical approach relies on the observation that the holomorphic Baxter equation (|2.11|) 
resembles a one-dimensional discrete Schrodinger equation. A specific feature of the model is that 
X— coordinates entering the Baxter equation takes complex values ( p.lOj ) and the Planck constant 
equals unity h = l.Q As a consequence, the quasiclassical limit corresponds to large values of the 
"energies" q2, ■ ■ ■ ,qN in Eq. (|1.10|) . 



2 



In quantum models, like the Toda chain, the Planck constant controls a shift of the argument of the 



Q— function in the l.h.s. of the Baxter equation (2.11) 
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To perform the large qr— limit in (|2.11|) , one introduces an arbitrary auxiliary parameter rj 
and rescales simultaneously the coordinates, x x/r], and the charges, qk QkV^- Defining 



tN{x) = V tN{x/v) = + q2X^-' + ... + Qm , (3.1) 

with On = QuV"' = 0{ri^) as 77 ^ 0, one finds that this transformation allows one to get rid of 
large parameters in the Baxter equation ( |2.11| ). At the same time, it sets the Planck constant as 



h = rj. Let us look for the solution to the holomorphic Baxter equation in the WKB form |Tg, 



Q{x/r]) = exp ( ^ / dx S'{x) ] , S{x) = So{x) + 7]S^{x) + 0{rf) , (3.2) 



where S'{x) = dS{x)/dx and Xq is an arbitrary reference point. Its substitution into (|2.11|) leads 
to the following relations 



2cosh^'(x) = ^((a;) = -(lnsinh^'(a;))' + ^. (3.3) 

x" 2 X 



One can systematically improve the WKB expansion ( p.2| ) and express subleading corrections to 
S{x) in terms of the leading term Sq{x). Similar expressions can be obtained for solutions to the 
antiholomorphic Baxter equation Q{x). Then, a general WKB expression for Q{x, x) is given by 
a bilinear combination of Q{x) and Q{x). 

3.1. Properties of the WKB expansion 

Introducing notation for 

= S'q{x) , y{x) = 2x^ sinhp^ , (3.4) 
one rewrites the first relation in ( |3.3| ) as 

y^{x) =t^{x) -Ax^^ . (3.5) 

We notice that, up to rescaling of parameters, y{x) coincides with the hyperelliptic curve F^v, 
Eq. ( 1.10| ). The coincidence is not accidental of course. The leading term of the WKB expan- 



sion, So{x), satisfies the classical Hamilton- Jacobi equations in the separated coordinates. Its 
derivative, Px = Sq{x), defines a (complex- valued) holomorphic component of the momentum in 
the separated variables. As such, it belongs to the spectral curve of the classical model ( |2.7| ) for 
w{x) = x^ exp^px). 

The leading term of the WKB expansion ( p.2|) can be calculated as 

So{x) = / dxpx= / —-^\NtN{x)-xtj^{x)\+xpx . (3.6) 

Jxo Jxo y{^) ^0 

Solving ( ^75|) we find that y{x) and, as a consequence S'q{x), are double- valued functions on 
the complex x— plane. To specify two branches of Sq{x), one makes cuts on the x— plane in an 
arbitrary way between the 2{N — 1) branching points aj. The latter are defined as y{(Jj) = 0, or 
equivalently 

?^(a,) - 4af = (g^af-^ + ... + g^)(4af + g^^"' + •.• + g}v) = . (3.7) 
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According to their definition, the branching points correspond to the special points on the phase 
space of the classical system, in which the holomorphic component of the momentum (in the 
separated coordinates) takes the values Px = and = ±iiT. Two different solutions to ( |3.5| ) 
give rise to two branches S'o,+ (x) and S'o,_(a;) which are continuous functions of complex x except 
across the cuts. These functions are transformed one into another as x encircles the branching 
point CTj in the anticlockwise direction 

S'oA^) -S'o,^{x). (3.8) 

Their asymptotic behaviour at infinity can be found from ( p.6|) and ( |3.5| ) as 

^1/2 

5^,±(a;)~±^~±^r^(/^-l/2), (3.9) 



as X oo. Here, in the last relation, we replaced % = q2'rf by its expression, Eq. (|1.5|), and took 
the limit 77 — with \ri{h — 1/2) | = \ri{ivh + n/i/2)|=fixed. Notice that the integration contour 
in ( p. 61 ) does not cross the cuts. 

It becomes convenient to combine the two branches S'q^{x) and define Sq{x) as a single- valued 
function on the Riemann surface Fat obtained by gluing together two copies of the x— plane 
along the (A^ — 1)— cuts [o"2j; o"2j+i] running between the branching points. By the definition, 
S'q{x) = S'q^{x) on one plane (upper sheet) and S'q{x) = S'q_{x) on another one (lower sheet). 
Then, it follows from Eq. ( p.6| ) that dx Sq{x) is a well-defined, meromorphic differential on Fjv of 
the third kind (the dipole differential). It has a pair of poles located above the point a; = cxd on 
the upper and lower sheets, P+ and P^, respectively. Here we used the standard notation for 
the points on the Riemann surface, = {x, ±). 

Substituting (|3.(j|) into (|3.3|), we calculate the first nonleading WKB correction as 



i f y{x)\' iNs ' ^ 



2 V""" 2x^ J ' X 



4^x-a,- X \ 2 J 

7 = 1 -J N ' 



^U^) = o +-- = 7 >. — - + — • (3-10) 



In distinction with the leading 6*0— term, the function S[{x) is well-defined on the complex 
x— plane. Therefore, it takes the same value on the both sheets of the Riemann surface 
and its asymptotic behaviour for a; ^ 00 is given by 

S[{x) -(ns-IY (3.11) 



x\ 2 



Combining together (^]6|) and (|3.10| ) we find that the two different branches S'q. 



X , or 



equivalently two different sheets of the Riemann surface Fat, give rise to two independent WKB 
solutions to the holomorphic Baxter equation (|2.11| ) 



Q±{x/r])=exp(- [ dxS'^ix)), S'^ = S'^^^^x) + r] S[ix) + 0{r]^) . (3.12) 



Their asymptotics at infinity can be obtained from (|3.9| ) and (|3.11|) as 



Q+{x/r]) ~ x'-'^-''^ , Q-{x/v) - x'^-''^ , (3.13) 
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as X — > oo. Going over through similar analysis of the antiholomorphic Baxter equation, one 
arrives at the WKB expressions for Q^{x/ri). They can be obtained from (|3.12|) by replacing 



holomorphic variables by their counterparts in the antiholomorphic sector. In this way, one gets 
from (|3l3| ) 

Q4x/r]) ~ x^-^-^^" , g_(x/r/) ~ x^-^^" , (3.14) 
where x = x*, s = 1 — s* and h = 1 — h*. 



3.2. Quantization conditions 



Let us construct the quasiclassical solution to the Baxter equation ( |2.11D as a bilinear combination 
of the chiral solutions Q±{x/ri) and Q^{x/rj) 

Q{x/ri,x/ri) = c+Q+{x/r])Qj^{x/r]) + c_ Q-{x/r])Q_{x/r]) . (3.15) 

Using (|3.13|) and (|3.14|) one verifies that the wave function defined in this way has correct asymp- 
totic behaviour at infinity, Eq. ( [2.121 ). The cross-terms Q±Q:^ do not enter ( |3.15| ) since they do not 
verify (|2.12| ). The functions Q±{x/ri) and Q^{x/r]) depend on the reference point Xq, Eq. (|3.12|) , 
whereas Q{x /rj,x /rj) should not depend on the choice of Xq. This fixes the Xq— dependence of 
the coefficients c±{xq) as (xq = Xq) 



c±[Xq) = c±(xojexp - 
1 ^ 



dx S'Ax) + 



dx S_i 



X 



xo 



(3.16) 



By the construction, the WKB formula (|3.15|) is valid for small t] and x ^ rf* . It describes the 
wave function, Q(x, x), for large values of the separated coordinates, x ~ x ~ l/?7, or equivalently 
n u l/rj m. the parameterization ( p.lO|) . In this region, one can ignore the fact that n 
takes strictly integer values and treat the separated coordinates x and x as continuous complex, 
mutually conjugated variables, x = x*. We recall that in one- dimensional lattice models, the 
Toda chain [ITB], |TB[ and the S'L(2,M) magnet [ITU], EDI, pT[, the classical motion in the separated 
coordinates is restricted to finite intervals on the real x— axis. The WKB wave function is a 
continuous, single-valued function of real x, oscillating inside these intervals and vanishing at 
infinity. Going back to the noncompact SL{2, C) magnet, one finds that, in distinction with 
the models mentioned above, the classical motion in the separated coordinates occurs on the 
whole two-dimensional x— plane. This suggests that the WKB wave function Q{x,x = x*) has 
to be well-defined on the complex x— plane. In particular, contrary to the chiral solutions to the 
Baxter equation ( |3.8| ), it should have a trivial monodromy around the branching points aj. In 
other words, Q{x,x*) has to be a single- valued function on the complex x— plane rather than on 
the Riemann surface F^v. The former condition is much stronger than the latter one and, as we 
will show below, it leads to the WKB quantization conditions for the integrals of motion q. 

Let us examine the monodromy of the chiral solutions to the Baxter equation, Q±{x/rj), 
around the branching points aj, Eq. ( |3.?1 ). Encircling the branching point aj on the x— complex 



plane, one finds that the leading WKB term S'q(x) is transformed according to (|3.8| ) while the 
subleading term S[{x) stays invariant. Let us explore a freedom in choosing the reference point 
Xo in ( p.l2| ) and put Xq = ctj in order to ensure that S±{aj) = 0. Then, it follows from (|3.8| ) and 
( p.l2| ) that the WKB solutions Q±{x) defined in this way are transformed one into another as x 
encircles aj on the complex plane 



Q±{x/t]) Q^{x/r]), (forxo 



(3.17) 
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Similar relations hold for the antiholomorphic solutions Q^{x/ri) at Xo = cxj = a*. We find from 
( p.l5|) that Q{x/ri,x/ri) stays invariant under this transformation provided that 



1,2, 



,2(iV-l)) 



(3.18) 



These conditions ensure that the quasiclassical wave function (|3.15|) is a single- valued function on 
the complex x— plane. We recall that the branching points are defined as solutions to Eq. ( p.7| ). 

For different Xq, the coefficients c±{xq) are related to each other according to (|3.16|) . Therefore, 
choosing Xq = aj and = we obtain from Eqs. ( |3.16| ) and ( |3.18D 



^-expf^Re 



dx [S'q _ [x] — S'q 



\x 



0{ri) 



(3.19) 



where j = 1, ...,2(A^ — 1) and the integration contour does not cross the cuts on the complex 



a:— plane. In arriving at ( p.l9|) , we applied (|3.12|) and took into account that S'q _|_(x) = [Sq aJ^x^* . 



Notice that the exponent in the r.h.s. of ( p. 191) does not receive the (9(?7°)— correction, since S'^ix) 
is single-valued on the x— plane. Since the l.h.s. of ( |3.19|) does not depend on j, one gets from 
( p.l9|) the set of consistency conditions 



exp 



^Re 



dx{S',^_{x)-S',^^[x)\^0{ri) 



(3.20) 



The quantization conditions ( p.l9| ) and ( |3.2CI| ) can be expressed in a concise form in terms of 
the contour integrals on the Riemann surface Fat, Eq. (|1.10|) . We recall that the two branches 
5'q _|_(x) define the dipole differential dx Sq{x) on Fat. This allows one to rewrite ( |3.19D as 



MO) 
c_(0) 



exp 




dx S'q{x) + 0{r]) 



-2ie 



(3.21) 



with some (real) constant introduced for later convenience. Here integration goes over an 
arbitrary path on Fjv, which starts at the point Pq located above x = on the lower sheet and 
ends at the point Pq*" above x = on the upper sheet. 

Let us define the canonical basis of oriented cycles on Fat as shown in Figure |I[ The contour 
integral in Eq. ( |3.21| ) can be uniquely decomposed over this basis as 



p+ 



N-2 



N-2 



(3.22) 



with kj and irij arbitrary integer. Here, the path 7p-p+ goes from the point Pg lower 
sheet to the point Pq on the upper sheet and it does not cross the canonical cycles (see Figure |l]). 
The l.h.s. of ( |3.21| ) should not depend on the choice of the integration path, or equivalently on 
integers kj and rrij in Eq. ( p.22|) . This leads to 



Re 



dxS'Q{x)/r] = 7r^2fe-i 



Re (b (ixS'o(x)/?7 = 7r£2fc 
hk 



(3.23) 



12 



Figure 1: The canonical basis of oriented a— and /?— cycles on the Riemann surface Ftv- The 
dotted line represents the part of the /3— cycles on the lower sheet. The cross denotes a projection 
of the points onto the complex plane. The path 7p-p+ goes from the point Pq~ on the lower 
sheet to the point on the upper sheet. 



with 5*0 (x) defined in (|3.6| ), A; = 1, A^ — 2 and I = . . . , i2N-\) being integer. These relations 
establish the WKB quantization conditions for the integrals of motion of the model. In addition, 



one finds from (|3.22|) and ( p. 21 



Re 



dx S'Q{x)/rj = 9 



(3.24) 



The following comments are in order. 

Eqs. ( ^21 and ( ^21 are valid up to corrections suppressed by a second power of 77. One 
can show that this feature is rather general and the WKB expansion in the l.h.s. of ( |3.23|) and 
( p.24|) goes over even powers of i]. 

We recall that, by the definition, rj is a small auxiliary parameter which was introduced in 
( |3.1|) in order to formulate the WKB expansion. The spectrum of the model should not depend 
on rj. Indeed, one verifies using ( p.6[ ) and (^.11), that the l.h.s. of ( p3.23D and ( |3.24| ), as well as 
Qn = Qn/f]'^-, stay invariant under the scaling transformation x — > Ax, rj ^ Xr] and g„ ^n^". 
Therefore we may choose rj in Eqs. ( |3.23| ) and ( ^.24 ) to our best convenience, say rj = 1, but keep 
in mind, of course, that the WKB approximation is valid for large values of the charges g„. In 
this way, one arrives at the quantization conditions ( |1.13| ). Their solutions have the form ( |1.4| ) 
and are parameterized by the set of integers £ = (£1, ^2(Ar-2))- 

As we will show in Section 4, the phase G in the r.h.s. of (|3.24|) is closely related to the 
quasimomentum On corresponding to the eigenstate (|3.15|) . Evaluating the contour integral in 
the l.h.s. of ( p.24|) and replacing the charges g„ by their quantized values, Eq. (|1.4| ), one can 
obtain from (p.24|) the dependence of the quasimomentum on integers £. 



A novel feature of the quantization conditions ( |3.23| ) compared to the conventional WKB 
approach |]T^ ^ is that they involve both a— and /?— periods on the Riemann surface T^. We 
remind that the Hamiltonian of the SL{2, C) magnet is given by the sum of two one-dimensional 
mutually commuting Hamiltonians "living" on the z— and ^— lines. If we z and z have been 
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real coordinates, each of these Hamiltonian would have defined quantum SL{2,M.) Heisenberg 
magnet. The WKB quantization conditions for this magnet look as follows fl^, ^ 



SL(2, M) : (p dx SQ{x)/ri = 27r(4 + 1/2) , (3.25) 



with the "action" differential, dx S'q{x), and the Riemann surface, T^, the same as in the SL{2, C) 
case, Eqs. (|3.6| ) and ( |1.10| ), respectively. A crucial difference between the SL{2, C) and SL{2,M.) 
magnets is that the integrals of motion, and the separated variables, x and px, take real values 
in the latter case. As a consequence, the classical motion in the separated SL{2, M) coordinates 
is restricted to the intervals on the real x— axis on which tj^ix) — 4x^^ > (see Eq. ( |3.5| )). 



The integration in ( |3.25| ) goes along the cycles (real slices of F^r) encircling these intervals on 
the complex x— plane. In the SL{2, C) case, the classical trajectories wrap arbitrarily around F^r 
and, as a consequence, the quantization conditions ( |3.23| ) involve also /?— cycles which correspond, 
from point of view of the SL{2,M.) magnet, to classically forbidden zones. 



We shall solve the quantization conditions ( 3.23|) in Section 5. 



4. Quasiclassical spectrum 

The WKB analysis performed in the previous Section allowed us to formulate the quantization 
conditions for the integrals of motion, Eqs. (|3.23|) , and construct the WKB expression for the 
wave function in the separated coordinates, Eq. ( |3.15|) . Let us extend our analysis and evaluate 



the physical observables of the model - the energy, E^, and the quasimomentum, 9]\f. 

Our starting point will be the expressions for En and 9n obtained in Ref. |10] within the 
method of the Baxter Q— operator (see Eqs. ( [4.6|) and ([4.7|) below). They are formulated in terms 
of the eigenvalue of the Baxter Q— operator, Q{u, u), which is a function of two complex spectral 
parameters, u and u, such that i{u — u) = n with n arbitrary integer. The wave function in the 
SoV representation, Q{x,x) (see Eq. (|2.9D), coincides with this function for m = z/ — m/2 and 
u = V — in/2 with v real. In this way, Q{u, u) can be considered as an analytical continuation of 
the wave function from the real axis to the whole complex z/— plane. 

The energy spectrum of the quantum SL{2, C) magnet is related to the behaviour of the 
function Q{u,u) around two special points on the complex i/— plane corresponding to u = ±is 
and u = ±is with (s, s) being a single-particle 5'L(2,C) spin, Eq. (|1.1| ). We notice that this 



behaviour can not be deduced from the obtained WKB expression for the Q— function (|3.15|) 
since the latter is valid only for large u and u. In this Section, we shall construct an asymptotic 
expression for Q{u,u), valid for large charges g„ and u,u = fixed, and use it to calculate the 
quasiclassical energy spectrum. 

4.1. Baxter Q— blocks 

To begin with, we summarize, following [l^, the main properties of the eigenvalues of the Baxter 
Q— operator, Q{u,u). The function Q{u,u) can be decomposed into a bilinear combination of 
chiral blocks 

Q(u, u) = e-^' Qoiu)Qoiu) - e'' Qiiu)Q,iu) , (4.1) 



14 



with 6 arbitrary complex. The blocks Qo{u) and Qi{u) satisfy the chiral Baxter equation (|2.11|) 
for arbitrary complex u and fulfil the Wronskian condition 

r^du - s) 

Qoiu + t)Q,{u) - Qoiu)Qiiu + t)= ). . (4.2) 

1 [tu + s) 

Similar blocks in the antiholomorphic sector, Qq{u) and Q-^^{u), are defined as 

The Q— blocks are meromorphic functions on the complex plane. Their analytical properties can 
be summarized as 

Qoiu) = r^-i(l - s + zu)f{u) , Q,iu) = r^-i(l - s - iu)f{u) , (4.4) 

where f{u) and f{u) are entire functions. As was already mentioned, at u = x and u = x* with 
X given by (|2.10|) , Q{u,u) defines the wave function in the separated coordinates, Eq. ( ^^ . The 
constant 6 in (|4.1|) is fixed by the requirement that Qlu — in/2, v + i?7,/2) should not have poles 
for real v. The same condition can be expressed as |jlO|| 



Q(z(l-s) + e,-zs + e) = C(e°), Q(-zs + e, z(l - s) + e) = 0(e°) , (4.5) 

for e — s> and s = 1 — s*. 

The energy and the quasimomentum of the model are expressed in terms of the Q— blocks as 

En = -2lm{\nQoiis)y + 2lm{\nQoi~is))' + En , (4.6) 
On = ,ln^^lMSi!!L. (4.7) 

where En = 2NRe [ip{2s) + ^{2 - 2s) - 2^{l)\. 

Thus, the problem of calculating the energy spectrum of the model is reduced to finding the 
blocks Qq{u) and Qq{u). Their exact expressions were obtained in ||10|. In this Section we shall 
obtain asymptotic expressions for the blocks Qq{u) and Qq{u) valid for large charges g„ and fixed 
spectral parameters u and u. 

4.2. Asymptotic solutions to the Baxter equation 



Let us rewrite the Baxter equation (|2.11|) for the block Qq{u) as 



{u + is)''q{u) + {u- isf , ^ = tN{u) , (4.8) 

q[u — i) 

where the notation was introduced for q{u) = Qq{u + i)/Qo{u). We notice that for large charges 
Qn and fixed u the polynomial t^iu), defined in ( p..lO| ), takes large values \tN{u)\ ^ 1. This 
suggests to expand the solutions to (^l8| ) in inverse powers of tN{u). Assuming that one of the 
terms in the l.h.s. of (14.81) is much smaller than another one, we obtain two solutions 



q^{u) = j^^^ + ..., q.{u-^) = ^-^ + .... (4.9) 
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Here, ellipses denote subleading corrections controlled by the parameter 

(m±^s)^(m±?(1-s))^ 



tN{u)tN{u ± i) 



< 1. 



(4.10) 



Eq. ( |4.9| ) gives rise to two linear independent asymptotic solutions for the block Qo{u). The 
general expression for Qo{u) with the prescribed analytical properties ( [4 .41) takes the form 



r 



N-li 



1 — s + iu] 



T{s-iu) ■ T^{s + iu) 

where if±{u) are entire functions satisfying the functional equations 



(4.11) 



(4.12) 



It is important to realize that the particular form of the ratio of the F— functions in the r.h.s. of 
( [4.11|) is uniquely fixed by analytical properties of the block Qoiu), Eq. ( ^^ . If any F— function 
in (^4.11 ) was substituted as T{x) — > 1/F(1 — x), one would obtain another solution to the Baxter 
equation (^4.8|) (up to corrections suppressed by the factor (|4.1CI| )) but its pole structure would be 
different from ( [4.4| ). 

The Q— block in the antiholomorphic sector is given by similar expression 



Qo (^) 



F 



N-li 



1 — s — iu) 



V9+(m) + 



1 



F(s + m) ""'"^ ' F^(s-m) 
with '^±{u) entire functions satisfying the relations 



(4.13) 



(4.14) 



Here tj\f{u) is given by ( |1.10[ ) with the charges g„ replaced by their antiholomorphic counterparts 
Qn = Qn^ so t^^^ tN{u) =tiy{u*)*. Substituting ( [4.1 1|) and ( [4.13| ) into (^) we obtain 



(as). . sin(7r(s + m))^ 
Qi (u) = F [1 - s + iu){ip+{u )) + 

TT 



sin(7r(s — iu)) 



T^{s + iu) 

F-(l - s - ^u){vA^*)r + TW7^-^(V-my 



TT F" (s — iu^ 

To solve ( [4. 12|) and ( |4.14| ) one factorizes the polynomial t]\f{u), Eq. ( [1.10| ), as 



(4.15) 



N 



tN{u) = 2u^ + q2U^-^ + ... + gAT = 2 JJ(m - Afc) . 



k=l 



For arbitrary, large g„ ~ 1/?]" its roots take, in general, large complex values, 
satisfy the sum rules 



^Afc = 0, AfcAn = g2/2 . 



Ha. 



-l)V/2 



(4.16) 
l/?7 and 

(4.17) 



k>n 
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Substituting ( [4.16|) into (|4.12| ) one can write a particular solution for (p+{u) in the form 



N 



(^V""™) {u) ~ e±™ 2-^" Yl r(zAfc - tu) . (4.18) 



k=l 



As before, substituting r(a;) l/r(l— x) in the r.h.s. of ( [4.1^ ) one can get yet another solution to 
( [4.12|) . To fix this ambiguity we require that ^p+{u) has to be an entire function of u in the region 



u ~ 1]^. Therefore, the product of the F— function in the r.h.s. of ( |4.18|) should not generate poles 
for u ~ 1]^. Decomposing into real and imaginary parts, T{iXk — iu) = r{i Re — Im Ajt — iu), 
one finds that in spite of the fact that A^ ~ l/?7, the roots with ImAfc ~ l/rj and ReA^ ~ rj^ 
generate the sequence of poles located at iu = i Re A^ — [Im A^] + ~ 77° with n nonnegative 
integer and [...] denoting the entire part. 

This suggests to separate all roots in ([4.18|) into two sets according to Im A^ > and Im A^ < 



and look for the solutions to ( [4.12| ) in the form 

ip+{u) = a+{u)if+{u), = a_(M)^_(M) , (4.19) 

where the notation was introduced for the basis functions 



r(l — iXu + iu) 



^-H = 2- n r(-.A. + .«) n rn-.n,-.„v (4.20) 



r(l + iXk — iu) 

ImAfc>0 ImAfe<0 V ' « ^ 



In these expressions, the roots with Re A^ ~ 77° do not generate spurious poles. In similar manner, 
the solutions to the antiholomorphic relations in ( |4.14| ) are given by 



If^u) = a+iu) (^+(r))* , If^u) = a„(n) (^4«*))* ■ (4-21) 



Substituting ([4.19|) and ([4.21|) into ([4.12|) and ([4.14|) , respectively, we find that the functions 



a±{u) and a±{u) are entire (anti)periodic functions satisfying the functional relations 

a±{u±i) = T{-l)^-a±{u), a±(M ^ i) = T(-l)^-a±(M) . (4.22) 



Here, N_ denotes the number of roots ([4.17|) with Im A^ negative 



= 1 , N+= 1 , N+ + N_ = N. (4.23) 

ImAfc<0 ImAfc>0 

Inserting ( ^.19| ) and ( ^.21| ) into ( ^.11| ) and ( [4.13| ), one obtains asymptotic expressions for the 
Baxter blocks Qo{u) and Qo{u). They depend however on four yet undefined functions a±{u) 
and a-j-{u). Additional constraints on these functions are derived in Appendix ^ (see Eq. ( |A.14] )). 
They come from two different sources. First, one has to ensure that the obtained expressions for 
the blocks verify the Wronskian condition ( [4.2|) . Second, for u = v — in/2 and u = u* (with u real 
and n integer) the function Q{u, u) coincides with the wave function in the separated coordinates. 
Substituting the asymptotic expressions for the blocks into ( |4.1|) and continuing the resulting 
expression to the region of large m, one should be able to match it into analogous expression for 
the WKB wave function ( p.l5|) . The matching procedure is performed in Appendix 0. 
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4.3. Energy spectrum 

According to ( [4 .61) and (|4.7| ), the energy and quasimomentum are related to the behaviour of 
the blocks Qo{u) and Qo{u) around the points u = ±is and u = ±is. Substituting the ob- 
tained asymptotic expressions for the blocks, Eqs. ( [4 .111) and ( [4.13| ), into the expression for the 
quasimomentum, Eq. ( |4.7| ), one finds after some calculation (see Appendix ^ for detail) 



e 



N 



-20 = -2 Re 



X] 



-2 Re 



p+ 



dx SLi 



(mod 27r) 



(4.24) 



Here in the second relation we replaced the phase 6 by its expression ( p.24| ) and put rj = 1. The 
contour integral in the r.h.s. of (|4.24| ) depends on the integrals of motion g„. Since the possible 
values of the quasimomentum are given by 6'7v = 27ii/N with i being integer, one should expect 
that for Qn satisfying the quantization conditions ( p.23| ) the integral takes the same quantized 
values. Indeed, we demonstrate this property in Appendix by an explicit calculation of ( 4.24 ). 

The calculation of the energy ( |4.6| ) goes along the same lines (see Appendix 0). It leads to 
the following asymptotic expression for the energy 



E^^'^ = 41n2 + 2Re ^ [^/'(l - s - iA^) + ^(s - iA^) - 2^(1) 



ImAi.>0 



+ 2 Re J2 [^C^ - s + iXk) + + tXk) - 2^{1) 

Im Afe<0 



(4.25) 



Here, the sum goes over the (complex) roots of the polynomial tN{u), Eqs. ( [4. 161 ) and ( |4.17] ). By 
the definition, their total number equals the number of particles, A^, and their values depend on 
the integrals of motion g„. Since the latter are quantized according to ( |1.4| ), the relation ( |4.25| ) 
establishes the dependence of the energy on the total set of quantum numbers, Eq. ( |1.4| ). 

The obtained expression for the energy ( |4.25| ) is symmetric under s — 1 — s. This is in 
agreement with the fact that the SL{2, C) representations of the spin s and 1 — s are unitary 
equivalent and, as a consequence, the corresponding spin magnets should have the same energy 
spectrum. 

According to (|4.25|) , the roots with ImA^ > and ImA^ < provide different contribution 
to the energy. To elucidate this property let us examine ( |4.25| ) for the SL(2, C) spin s = 



N 



E 



(as) 
N 



41n2 + 2Re^ U(l+iReAfc + |ImAfc|)+V^(iReAfe + |ImAfc|)-2V^(l) , (4.26) 



fc=0 



where A^ are roots of the polynomial t^lu), Eq. ( [4.16|) . We recall that the Euler ■?/'— function has 
poles at nonpositive integer values of its argument. In the r.h.s. of ( |4.26| ) these poles are never 
approached for arbitrary charges g„ provided that qn 7^ 0. 

Let us compare the asymptotic expression for the energy, Eq. ( |4.26| ), with the results of the 
1^, |lOl- For the sake of simplicity, we choose the total Lorentz spin of the 
0, so that the total SL{2,C) spin ( |1.5| ) equals h = 1/2 + iu^, and examine 

n(as) I 



exact calculation 
model to be Uh = 

the dependence of the energy E^^'\s=o on real Uh along a few lowest lying trajectories at = 3 
and A^ = 4. Applying ( [4.26| ) and ( |4.16|) , we substitute q2 = 1/4 + and replace the charges 
qs, . . . ,qN by their exact values found in Ref. R M. Comparison with the exact expression for 
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Figure 2: Dependence of the energy on the total spin h = 1/2 + iuh for the three lowest 

eigenstates at = 3 (left panel) and = 4 (right panel). The solid lines denote the exact 
energy from Refs. ||^, |l^, the dotted lines describe the asymptotic expression ( |4.26| ) calculated 
for the exact charges gs and 54. 



the energy is shown in Figure 0. We recall that Eqs. ( f4.25| ) and ( [4.26D were obtained to the 
leading order of the WKB expansion under assumption that the "energies" q2, ■ ■ ■ ,qN are large. 
Therefore, it is not surprising that the quasiclassical expression for the ground state trajectory 
agrees with the exact energy only for > 1. At the same time, for the excited trajectories the 
agreement is rather remarkable (especially at A^ = 4) even for smaller u^. 

In Section 3.2 we already drew the analogy between the quantization conditions for the 
SL{2,C) and SL{2,'R) magnets. It can be further extended to the asymptotic expressions for 
the energy. For the SL{2,M.) magnet the corresponding expression looks like |T9|, ^ ^ 



SL{2, 



E 



(as) 
N 



N 

2 In 2 + Re ^ [^(1 + iA^) + V^(iAfe) - 2^(1) 



(4.27) 



k=0 



Here, in distinction with the SL{2, C) case, the roots take strictly real values. Similarity 
between Eqs. ( |4.26| ) and ( |4.27|) suggests that there should exist a relation between the energy 
spectrum of the SL{2, C) and SL{2, M) magnets. Indeed, it can be shown that the latter can be 
obtained from the former by analytical continuation in the total spin of the system, Eq. (pTS]), 
from real to pure imaginary u^. 



5. Solving the quantization conditions 

In previous Section we derived the WKB quantization conditions for the integrals of motions 
of the SL{2,C) magnet, Eq. ( p.23| ), and obtained the expressions for the energy spectrum, 
Eqs. ( [4.25|) and ( |4.24D . In this Section, we shall solve the quantization conditions (|3.23| ) and 
reconstruct a fine structure of the spectrum. 

The general solutions to (|3.23| ) have the form ( [L.4|) . The spectrum of quantized charges qn is 
parameterized by the set of integers £ = (£1, . . . ,i2N-4) entering the r.h.s. of ( p.23| ), as well as 
by continuous real and integer nh defining the total SL{2, C) spin of the magnet, Eq. ( |1.5|) . 
We recall that the quantization conditions ( p. 231 ) involve the periods of the "action" differential 
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over the canonical basis of oriented cycles, a = (ai, . . . , aN-2) and f3 = {Pi, ... , Pn-2), on the 
Riemann surface Fat of genus g = N — 2. The definition of this basis on r^v is ambiguous ||2^ 



a 



a 



a ■ a. + b ■ 13 , f3^f3' = c- OL + d- (3, 



(5.1) 



where a, 6, c and d are (A^ — 2) x [N — 2) matrices with integer entries such that 

,M 1 



a b 
c d 



det Z = 1 



-1 



1 

-1 



(5.2) 



so that Z e Sp(iV— 2, Z). Notice that Z G SL{2, Z) for = 3. The spectrum of the model should 
not depend on the choice of the basis. Indeed, the quantization conditions (|3.23|) stay invariant 
under the Sp(A'" — 2, Z) transformation ( |5.1| ) provided that the integers ^odd = (^i, • • • ,^2N-5) 
and £even = (^2, • • • , ^2N-i) are transformed in the same way 



-odd 



odd 



a ■ £, 



odd 



b-L 



£' 



(5.3) 



This relation establishes the correspondence between two different solutions to the quantization 
conditions labelled by the sets of integers £ and £'. 




Figure 3: Flow of the quantized values of the integrals of motion q 



...,qN) with Uh. 



As was shown in Ref. [|^, quantized values of the charges, qn{i'',nh,£), form the family of 
one-dimensional nonintersecting trajectories on the (A^ — 1)— dimensional space of q = (g2, Qn) 
(see Figure Each member of the family is labelled by integers Uh and £, while the "proper 
time" along the trajectory is defined by Uh. The quasi classical approach allows us to reconstruct 
the trajectories in two different hmits: (i) £ = large and Uh = fixed; (ii) £ = fixed and Uh = large. 
In the both cases, the integrals of motion take large values, so that the quasiclassical approach 
is applicable. In the first case, we shall find the points on the q— space, at which the trajectories 
pierce the hyperplane Uh = fixed, and demonstrate that their positions define a lattice-like 
structure on the moduli space (see Figure |^). In the second case, we shall describe the flow of 
these points with along the trajectories. 
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5.1. Lattice structure 



Before we proceed with evaluating the contour integrals entering (|3.23|) , let us rewrite the quan- 
tization conditions in a slightly different form. ^From the expression for the quasimomentum 
( |4.24|) one gets 

r 27r 
9n = 2Re dxSo{x) = —i (mod 27r) . (5.4) 

Here the integration contour 7(0"^) goes from the point Pq*" above x = on the upper sheet of 
Fjv to the branching point cjfc, encircles it and goes to the point Pq" '^^ ^^e lower sheet of r^v- 
The quasimomentum Oj^j describes the transformation properties of the eigenfunction under the 
cychc permutations of particles, ^1/(^2, • • • , ^n, ^1) = e'^'^ ^(^i; • • • ; ^tv-i, z^)- Its value satisfy 
exp(zA^07v) = 1, so that i is integer in the r.h.s. of ( |5.4|) , < £ < — 1. One concludes from 
that 

If , £ 

-Re/ dx S'(){x) = — + Uk , (5.5) 



with k = 1, ...,2{N — 1) and integer. Remarkably enough, the quantization conditions ( 3.23|) 



are equivalent to the system of equations (|5.5|) . To see this, one rewrites the a— and /?— periods 
of the "action" differential as (see Figure pJ) f = f , ^ — f ^ ^ and <f « = f / x — f / x. 
Substituting these relations into (|3.23|) , one finds that the quantization conditions ( p.23| ) can 



be expressed as a linear combination of integrals entering the l.h.s. of ( p.5| ). In this way, one 
establishes the correspondence between two sets of integers 

n2k+i - = hk-i , rii - = hk ■ (5.6) 

Eq. ( ^.5|) involves a rather complicated contour integral on the hyperelliptic curve F^v- Although 
it is straightforward to calculate it numerically for a given set of the integrals of motion q, this 
does not allow us to understand a general structure of solutions to ( ^.51) . 
To this end, let us examine the quantization conditions (|5.5|) in the limit 



q2 = 0(e°) , g„ = 0{e~^^^~'^/') , = 0{e-^/') (5.7) 

with e <^ 1 and n = 3, ...,N — 1. This hierarchy corresponds to large q^, . . . ,qM and fixed g2- 
The main advantage of is that the spectral curve Fat, Eq. ( |1.10D , simplifies significantly 
and integration in (p.23| ) can be performed analytically. Indeed, after the scaling transformation 
X — *• (gAr/4)^/^x, y q^y the spectral curve (|1.1CI| ) takes the form 



fS^^) : y' = (x^ + 1 +p^_2(x))(l +p;v-2(x)) + 0{e^) , (5.^ 



where 



N-l 



PN-2{X) = J2^n X^-" , «n = f (^) '^^^^ = 0{e) (5.9) 



n=2 



4 V 4 



Comparing ( |5.8|) with (|1.10|) and ( |3.7|) , we find that A^ branching points of the curve F^*^"* are 
located at the vertices of the A^— polygon, whereas the remaining N — 2 points are located far 
from the origin on the complex plane 
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where m = l,...,N and k = N + 1,..., 2{N - 2). 

It becomes straightforward to evaluate the hypereUiptic integral in (|5.5|) for the first set of 
the branching points in ( ^.10| ) (see Appendix ^ for detail). Combining together ( [B.13| ), ( [B.14| ) 
and ( p.5|) we find the system of (A^ — 2) equations for the integrals of motion 



(uN^l-mUN) ■ B ^) - {Um^lU^r ■ B Q, = TT ^ e-'^^^'^-^W^ , (5.11) 



with m = 2, N — 2. Here -B(a;, y) = r{x)r{y)/r[x + y) is the Euler beta-function, the moduli 
Um were defined in (|5.9|) for 2 < m < — 1 and un = {(In/^Y^^ ■ One also gets the following 
relation for the quasimomentum 



N 

(mod A^) . (5.12) 



m=l 



These relations were obtained in the small— e limit and they hold up to corrections ~ e^/"^. Notice 
that the sums in r.h.s. of Eqs. ( ^.11|) and (|5.12|) depend only on a subset of integers ni, . . . , n2N-i 



corresponding to the first set of the branching points in ( |5.1CI| ). Comparing the same small— e 
asymptotics of the the both sides of ( ^.11|) one finds that rik ~ e~^/^. Then, Eq. ( ^.61) leads to 



Replacing in ( |5.11| ) the moduli m„ by their explicit expressions (|5.9|), one obtains the system 



of (A^ — 2) equations for the integrals of motion q^, . . . ,qN- Since the second relation in ( ^.11 



IS 



invariant under m N—m, the number of independent relations reduces to (A^— 1) real equations 
and, therefore, the system (|5.11|) is undetermined. The remaining quantization conditions follow 



from the analysis of the integral ( |5.5| ) for the second set of the branching points in Eq. ( |5.10| ). A 
straightforward calculation shows that the corresponding n— integers, entering the r.h.s. of ( ^.5|) , 
scale in the limit ( p.7|) as Uk ~ e^^"^ and induce subleading WKB corrections. 

The quantized values of the "highest" charge qn can be calculated from the first relation in 
( ^.11|) . One finds after some algebra the following remarkable expression 



ni + 2/N) 



TT A^-2 

where the notation was introduced for the Fourier series 



cot{7i/N)\Qin)\-' + 0{\Qin)\-') 



(5.13) 



N 



Q(n) = ^n.e-(2'=-i)/^, (5.14) 



=1 



and = 1/4:+ (z/^ + inh/2Y - N{us + m^/2)2 according to (P3|). 

The sum in the r.h.s. of ( p.l4| ) is invariant under simultaneous shift of integers, — + a. 
This transformation changes the value of i in Eq. ( |5.12| ), but leaves invariant the quasimomentum 
9]\f = exp{2'iTii / N) . Therefore, q]^^ and 9^ depend on the differences — n^+i, which in their 
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turn can be expressed in terms of the integers with a help of (|5.6|) . To obtain the corresponding 
expressions for qj^'^ and On one chooses the "gauge" rii = and substitutes in ( |5.14D and (|5.12|) 



n2k 



-2k 



n2k+i — ^2k~i — ^: 



2k 



(5.15) 



with fc = l,...,A^ — 2. It follows from ( |5.13|) that, to the leading order of the WKB expansion, q^j^ 
does not depend on the total SL{2, C) spin h = {1 + nh)/2 + iuh- The /;.— dependence enters into 
( [5 .131 ) through the nonleading correction, which becomes smaller as one goes to higher excited 
states with larger n^. 



5.1.1. Special case: = 3 

For iV = 3 we find from ( 5.13 ) and ( |5.12 ) the quantized values of the charge qs 

r3(2/3) 



1/3 
Q3 



271 



-[rii - 2n2 + ris) + — (723 - rii) 



+ 0{e 



1/2 ^ 



(5.16) 



and the quasimomentum £ = —rii — n2 — n^, with ^1,2,3 ~ e Here, for simplicity we did not 
include the nonleading correction ~ Using ( |5.6| ) one rewrites these relations as[| 



(4,4) 
^3(4,4) 



r3(2/3) 



2n 
2n 



\{^l + ^2) + ^^{h-^2) 



(£1 + 4) (mod 27r) . 



(5.17) 



Thus, to the leading order of the WKB expansion, the quantized charges q^^'^{£i,£2) are located 
on the complex plane at the vertices of the lattice built from equilateral triangles. 



0.5 



is^ 



-0.5 



V '4 \¥r ''4—- 
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V]^ ^ 



-1 -0.5 0.5 

l/3i 



1 



1 /3 

Figure 4: Lattice structure at = 3. Crosses denote the exact values of gg at g2 = 1/4. Dotted 
lines intersect at the points defined in Eq. ( |5.17D . 



As one can see from Figure |^, Eq. (|5.17| ) is in a good agreement with the exact results. Notice 
however that the exact values of q^ do not approach the origin, so that a few lattice vertices 

•^To make the correspondence with Ref. one has to redefine integers as £1 + ^2 ^ and £1 — ^2 ~> ^2- 
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remain vacant. Since the corresponding gs are small, one should not expect the WKB approach 
to be applicable in this region. Indeed, one can verify that for | gal < A3, the first nonleading 
WKB correction to the periods 0.(^3) and a/? (53) becomes comparable with the leading order 
contribution. 



5.1.2. Special case: N = A 



For = 4 the spectrum of the magnet is parameterized by two quantum numbers and q^. 
From ( p.l3| ) one gets the charge q^ as 



94 



1/4 _ r^(3/4) 

4v^ 



1 i 
-^(^1 - 77-2 - ^3 + ^4) + -^{-ni - ^2 + n3 + n^) 



+ 0{e 



l/2^ 



(5.18) 



where, in general, '^'1,2,3,4 = C(e ^/^). The quasimomentum ( |5.12| ) is equal to 

64 = i = —{rii + n2 + n3 + n^) (mod 2tt) . 



To find the charge ^3, we apply the second relation in ( |5.11D at = 4, m 
definition of the moduli (h.9\) 



Im 



93 

1/2 

94 



-Til + n2 - ri3 + 714) + C(e 



3/2 ^ 



(5.19) 
2 and use the 

(5.20) 



As was already explained, the system ( 5.11 ) is undetermined and it does not fix the charge ^3 
completely. The additional relation on ^3 comes from the analysis of the branching points in 
( ^.10|) located far from the origin. 

The solutions to ( [5. 18] ) and ( |5.2(J| ) are parameterized by three integers ii = n^—n2-, ^2 = ni—n2 
and £4 = ni — n^. To reveal the properties of the spectrum it is more convenient, however, to 
introduce their linear combinations 



mi 
m2 



{rii - 77-2 - ris + n4)/2 = {rii + n^) + 



-rii -n2 + n3 + ni)/2 = {n^ + ^4) + - 



(5.21) 



where integer i defines the quasimomentum (|5.19|) . Notice that mi, 2 are integer for i = even and 
half-integer for i = odd. Choosing the gauge ni = 0, one rewrites (|5.18| ) and (|5.20|) as 



94 



/4^r^/mi m2| 1/2 
2v^ VV2 V2J ^ ' 



Im 



93 

1/2 

94 



2 ( mi - m2 - - 



0(e 



3/2N 



(5.22) 



where mi = (—£1 + 2^2 — ^4)/2, m2 = {^l — ^4)/2 and mi, 2 = 0{e ^/^). Let us examine these 
expressions in more detail. It is convenient to consider separately the iV = 4 eigenstates with 
53 = and gg ^ 0. 

For the eigenstates with gs = one finds from ( p.20| ) that ni+n^ = n2+n4. As a consequence, 
their quasimomentum, 64 = —ttI/2 = 7r(ni + n^), takes the values ^4 = 0, vr (mod 27r), while 
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mi = rii — n2 and 1712 = — n2 are strictly integer. Thus, the quantized values of for the 
eigenstates with ga = are described by the first relation in ( |5.22D with mi.2 integer. They form 
a square lattice on the complex gy**— plane whose vertices are specified by the pair of integers 
(777.1,777.2). For mi ± m2=even, i = (mod 4) and the quasimomentum equals 64^ = 0. For 
mi ± m2=odd, i = 2 (mod 4) and the quasimomentum equals 64 = n. 
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Figure 5: Lattice structure at = 4. Crosses denote the exact values of for different and 
g2 = 1/4. Dotted lines intersect at the points defined in Eq. ( |5.22D with mi_2 integer (left panel) 
and mi 2 half-integer (right panel). 



For gs 7^ the eigenstates can be separated into two groups according to their quasimo- 
mentum, 6*4 = 0, TT and 9^ = ±7i/2. Let us visualize the solutions to (|5.22|) as points on the 



three-dimensional ^— space with the coordinates ^1 = Reg]'''^, ^2 = Img]^'^ and ^3 = Im (gs/g^^)- 



O4 = 0, tt: One finds from ( [5.21| ) and ( |5.19| ) that i is even and mi 2 are integer. According 
to ( ^.22| ), g4'''* does not depend on i. Therefore, choosing 2(mi — m2) — £ = 0, ±2, ±4, 
one finds that the solutions to ( |5.22| ) define an infinite set of identical square lattices in the 
^— space. These lattices run parallel to the (^1, ^2)— plane, as shown in Figure ^on the left, 
and cross the ^3— axis at ^3 = 0, ±2, ±4,.... Exact results indicate that the degeneracy 
between lattices with different ^3 is lifted by nonleading WKB corrections to g3 and g4. 



64 = ±7r/2: One finds from ( |5.21|) and (|5.19|) that i is odd and mi 2 are half-integer. The 
solutions to ( |5.22| ) define an analogous lattice structure in the ^— space. The charges ql^^ 
form square lattices on the (^1, ^2)~plane, shown in Figure § on the right, which are dual 
to the similar lattice in the previous case and have the coordinates ^3 = ±1, ±3, ±5, .... 



As follows from Figure ||, Eq. ( |5.22|) is in a good agreement with the exact results of Refs. .Q 
Surprisingly enough, Eq. ( p.22| ) works throughout the whole spectrum including the ground state. 
The latter is located at g3 = and g4 given by ( ^.22 ) for mi = 2 and m2 = 0. Similar to the A^ = 3 
case, a few lattice sites remain unoccupied in the small g4 region, where the WKB approach is 
not applicable. 



^At the same time, Eq. ( 5.22 ) invalidates the claim of Refs. 
values. 



that the charge (74 may take only real 
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Figure 6: Lattice structure on the complex g^^— plane is defined by linear integer combinations 
of the vectors ei, . . . , ejv- 



Going over to higher A^, we find that the lattice structure becomes more complicated. The 
leading order expression for the "highest" charge, Eq. ( |5.13| ), can be written in the vector form 



i/AT _ J{l + 2/N)_ 



N 



Y.rikek + Oie'/''), (5.23) 



k=l 



where = e ^'^^'^^ ^^/^ define unit vectors on the complex plane as shown in Figure ^. The 
quantized values of g]^^ are given by linear combinations of these vectors with integer (positive 
and negative) weights. Since ^^=1 6fc = 0, the charge q]^^ depends on the differences — n^+i, 
or equivalently on the integers, Eq. (|5.6| ). 

5.2. Whit ham flow 

In the previous Section, we solved the quantization conditions ( |3.23| ) in the limit ( |5.7] ), which 
corresponds to large charges q3,...,qN and fixed q2, or equivalently the total SL{2,C) spin 
h = {l + nh)/2 + ii'h. Let us now determine the dependence of the charges ( |1.4| ) on the continuous 
parameter u^. 

To simplify analysis, we choose a single particle spin in (|1.5D as s = and consider the limit 



i^h ^ nh- One finds from ( |1.5| ) that q2 = 1/4 + z/^ takes large positive values in this limit.Q Let 



us explore an ambiguity in choosing the WKB parameter r] in Eqs. (|3.1|) and (|3.23| ) and fix it as 



12 



-'" = {l/A + ul)-''\ qn = qnq7'\ (5-24) 



with n = 2, iV, so that q^ = 1. In this Section, we shall solve the quantization conditions ( |3.23|) 
and determine the dependence of quantized g„ on r]. As an example, we present in Figure ^ the 
dependence q^ = qsii^h] h,^2) at N = 3 for different trajectories. We will show below that the 



dependence of the charges is governed by the Whitham equations [27 



^One can also consider the limit nh'3> Vh, so that q2 — 1/4 — n\/A is negative. Similar analysis allows one to 
find the n/j— dependence of the quantmn numbers 
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Figure 7: The dependence of % = q^/c^^'^ on the total spin h = 1/2 + iuh along different 
trajectories specified by integers {ii,i2)- The charge takes real values along the (2,0)— and 
(3, 0)— trajectories, and pure imaginary values along the (1, —1)— trajectory. Cross denotes a 
point with ^3 = ±i/^/27. 



5.2.1. Whitham equations 

Derivation of the Whitham equations is based on the following property of the "action" So{x) 



Eq. (Km, p8 



d_ 

dr] 



So(x) 



XO 



or] 



dx dt^^x 
y{x) dr] 



N 



\r- dqk 

^ dri 

k=3 ' 



dxx 



N-k 



y{x} 



(5.25) 



where t^ix) = 2x'^ + x^ ^ + X]fc=3 9k ^ and y(x) was defined in ( p.5|) . Differentiating the both 
sides of the first relation in (|3.23|) and taking into account ( |5.25D , one finds after some algebra 

^ f dx x'^~^ 



Re 



J2 dvQj 

J=3 



7CI 



2k-l 



(5.26) 



The second equation in ( ^.23| ) leads to a similar relation with — > (3^ and ^2k- The l.h.s. 

of ( |5.26| ) involves a (unnormalized) differential of the first kind on F^r. Its a— and /5— periods 
can be parameterized as 

dx x^'^ 



oik 



dx x^ ^ 
y{x) 



2n [U-^ 



jk 



2tt [U-U 



jk 



(5.27) 



where j, k = 1, N — 2. Here r = [Tjfc(g)] is the Riemann matrix for the hyperelliptic curve F^v, 
while the matrix U = [Ukj{q)] defines the normalized holomorphic differentials (see Eq. ( |C.1| ) in 
Appendix 0). Both matrices depend on the charges g„ and are independent from the flow param- 
eter rj. In addition, the Riemann r— matrix is symmetric and has positively definite imaginary 
part [0. 

To solve (|5.26| ) one considers linear combinations Xk = ^^Ji dr^qN+i-j [U~^{q)]jk- They 
satisfy the relations 

N-2 

2 Re Xfc = e2k-i , 2ReJ2 ^jk = hk , (5.28) 
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whose solution can be written in the matrix form as 



X( 



where the vectors 



(•^cvcn ' '^odd) ("^even "^i 



t 

odd 



■rt) 



(5.29) 



2Imr^™ ^""^ ^ '^^'^ "^'^ ^2Imr' 

(•^odd) are built from integers (^2fe-i), and Imr = (r — r"l')/(2i) is a 



positively definite symmetric matrix. Finally, one replaces X{q) in (|5.29|) by its definition, puts 
rj = q2^^'^ and finds the system of Whitham equations 



3/2dqn 
oq2 



(•^even -^odd ' 



1 



Imr( 



■U{q) 



(5.30) 



N+l-n 



with = 3, ...,N. Notice that the r.h.s. of ( 5.30 ) is g2~iiidependent and it depends only on the 
charges g^. 

The Whitham equations ( ^.30|) involve the matrices f/(^) and r(g) defined in (|5.27|) . Their 
explicit expressions depend on the choice of the canonical set of the oriented cycles on Fat. As 
was already mentioned, the a— and /3— periods are defined up to transformation (|5.1|) , which 

acts on U{q) and r(g) as 



U {a + bT)-'U, 



hT){a + hr) 



-1 



(5.31) 



It is easy to verify that the Whitham equations ( p.30| ) are invariant under ( |5.31|) provided that 
.^odd and £even are transformed according to (|5.3|) . 

The Whitham equations ( [5.30[ ) describe the dependence of the charges, g„ = g„(i^/j,£) on the 
total spin of the system z/^. They define the flow of the quantum numbers g„ with Vh along 
the trajectory labelled by the integers i (see Figure To solve (|5.30| ) one has to specify the 
initial conditions for g„ (with n = 3,...,N) at some reference q2- They are provided by the 
expressions for the charges g„, Eqs. ( |5.13| ) and ( p.23| ), obtained in the previous section. We recall 
that Eqs. (|5.13|) and ( ^.23|) were obtained in the region of the moduli space ( pl7|) corresponding 
to q2 = fixed. The Whitham equations (|5.30|) allow us to evolve the charges g„ to arbitrary large 
values of q2. 



5.2.2. Whitham flow at = 3 

In the rest of this section, we shall present a detailed analysis of the Whitham equations at 
N = 3. Generalization to higher N is straightforward and can be performed along the same 
lines. 



At = 3 it is convenient to introduce notation for the periods of the action differential (3.6) 



a((h) = — (t dx 



«d(^3) = (f dX' 
27r J a 



y{x) 2tt Jp y{x) 

Here integration goes over the a— and /5— cycles on the elliptic curve F3, Eq. (|1.10| ) 



F3: y\x) = {x + q3){4x^ + x + q3)=Al[{x-a,). (5.33) 

3=1 

The genus of the Riemann surface defined by F3 equals g = N — 2 = 1. The quantization 
conditions ( |3.23 ) for the charge qs = q^/c^^'^ look like 



Re 0(53) 



2q 



1/2 



Rea£)(g3 



2q, 



1/2 • 



(5.34) 
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The Whitham equations ( ^.301 ) take the following form at = 3 



3/2<9g3 
<l2 



dq2 4Imr(g3) 
with the functions [/(gs) and T^q^^) defined from (|5.27|) and ( |5.25|) as 



(5.35) 



1 



a'{q?. 



27r 



dx 
y{x) 



(5.36) 



where a'($3) = da{q^) / dq^ and similar for a'^. We recall that Imr(g3) > for arbitrary ^3. 

Performing integration in the r.h.s. of ( p. 32 ), one can evaluate a{q^) and ar){q^) in terms of the 
elliptic function of the first and the second kinds. The resulting expressions for a{q^) and 00(^3) 
are analytical functions on the complex gs— plane with two cuts. The cuts start at the values 
of gs, for which any two branching points of the curve merge, aj = ak, and the integrand in 
( p.32| ) develops a pole. This happens for gs — 00. The remaining singular points correspond to 
zeros of the discriminant of the polynomial in the r.h.s. of ( 5.33| ), 



16 n( 

j>k 



\2 

aj - ak) 



-g3^(l + 27g3 



(5.37) 



In this way, one finds another three singular points on the complex g3— plane 



53, 



smg 



27 



27 



(5.38) 



Thus, the two cuts run on the complex g3— plane between these three points and g3 = 00. As we 



will show below, the solutions to the quantization conditions ( |5.34|) can be obtained in a closed 
form at the vicinity of these points. Let us determine the asymptotic behaviour of the functions 
a(g3) and aD(g3) around the singular points, gs = 0,—i/\/TI and 00. The behaviour around 
gjj = i/^/27 can be found by making use of the symmetry of the curve ( p.33| ) under x — >■ —x and 

ga -ga- 

As the starting point, one has to specify the a— and /?— cycles on the curve r3 (see Figure |I|). 
It is convenient to choose them in such a way that the a— and /3— cycles shrink into a point for 
g3 —i/\/27 and gs — > 0, respectively. Obviously, this choice is not unique and one can use 
another definition of the cycles. The resulting expressions for the periods a(g3) and aoiqa) are 
related to each other by the SL{2, Z) transformation (|5.1|) and 

At g3 
and (74 = 
[0-2, 0^3] 



the branching points are located along the imaginary axis at (Xi = (T2 = 0, as = ?/2 
= —i/2. According to our definition of the cycles. Figure |l], the a— cycle encircles the cut 
whereas the /3— cycle shrinks into a point. Calculation of ( [5. 3 2] ) leads to 



az)(0) = 0, 



a(0) 



dx 



V4x2 + 1 



2 

JO 



i/2 



dx 



V4a;2 + 1 2 



(5.39) 



To obtain the behaviour of the functions ^(gs) and aniqs) in the vicinity of gs = 0, one examines 
their derivatives with respect to gs. According to (|5.25| ), they are given by the a— and /?— periods 
of the holomorphic differential dx/y{x) on Fs, which can be calculated by the standard methods. 
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The details of the calculations can be found in |[T^. In this way, one obtains the asymptotic 
behaviour of the periods around gs = as 



a ($3) 



IT 

m + 



-qs [\n{iq3) - 1] + ... 



(5.40) 



where ellipses denote sub leading (9 (g|)— terms. 

At qs = —2/ a/27 the branching points are located at = a2 = i/vT2, cri = i/V^ and 
(74 = The /?— cycle encircles (J2 and cxi, while the a— cycle shrinks into a point (see 

Figure |I|). The periods ( p. 32] ) are given by 



1 ri/VT2 

a{-i/V27) = , aD{-t/V27) = - / 

^ J^/^/27 



dx 



In 2 



TT 



(5.41) 



Expanding the periods in the vicinity of gs = — z/V 27 one finds 

0(93) 



3(1-93/9) + -, 

= — + 7r(i - ^s/^) [ 1^ - ^3/9) - c] + 



(5.42) 



Vv27, c = 1 + ln(27/2) and ellipses denote C((l - gg/g)^) -terms. 



where q = — 

For g3 ^ 00 the branching points (Xfc move away from the origin and the spectral curve (|5.33|) 
can be approximated as = (x + q3){4:X^ + $3). Three branching points are located at the 
vertices of equilateral triangle, cr^ = (^3/4)^/^ e*'^*^^'^'"^''/^ with k = 1,2,3, and the last point at 
(J4 = — ^3. The calculation of the periods ( |5.32| ) leads to 



^(93) 
aoiqs) 



IJ3-+00 



93 



1/3 



2tt 



313(2/3) 



93 



1 /3 

where ellipses denote subleading 0(^3" )— terms. Notice that ( |5.43| ) can be written as 0(^3) = 
(/3 — /2)/(27r) and 00(93) = {h — -^2)/(27r), where the integral Ik was defined in ( |B.12D for 
N = 3. We already encountered the same elliptic integral in Section 5.1, when we analyzed the 
quantization condition in another region of the q— space, Eq. ( ^.7|) . Matching (|5.24| ) into (|5.7|) 
we find that 

q, = Oie~'/') , q, = 0(e°) . (5.44) 

Thus, the two regions, Eqs. (IsTTI) and ( ^.24|) , overlap as gs — 00 and g2 = fixed. As a consequence, 
one expects that at large q^ the solutions to the quantization conditions ( |5.34| ) should match the 
expressions for the quantum numbers gs obtained in the previous Section, Eq. ( |5.23|) . Indeed, 
substitution of ( |5.43| ) into ( p. 341 ) yields the expression 



1/3 



27T 



313(2/3) 
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(5.43) 
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1/3 _ -i/2r3(2/3) 



92 



271 



+ 0{ql^') 



(5.45) 
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1/3 -1/2 
% <i2 ■ 



1 /3 

which coincides with (|5.171 ) since ^3 

Let us substitute the obtained expressions for the a— and a^i— periods into the quantization 
conditions ( |5.34|) and derive the WKB expression for the charge ^3. We remind that Eqs. ( |5.4CI|) 
and ( p.42| ) hold in the vicinity of ^3 = and % = —i/ v^27, respectively. 
For 1^3 1 <^ 1 one finds from ( |5.4CI| ) and ( [5 .341 ) 



' 2 ^2 5'2 5 



Re[g3(ln(2g3)-1)] 



TT 



g^iy2 1 



(5.46) 



where 



Eq. (|5.46|) defines the scaling behaviour of quantized gs in the region ^3 ~ q2 

2) by another pair of integers 



for large q2- Notice that in comparison with (|5.34|) we replaced (^1 
(£'^,£'2). This was done in order to distinguish the integers entering the r.h.s. of ( p.45| ) and 
( |5.46|) . As was already mentioned, the periods a{q^) and ao^qz) depend on the definition of the 
a— and /3— cycles on the Riemann surface r3. These cycles encircle the branching points cxfc (see 
Figure |I|), which are moved on the complex plane as q^ varies. The two pairs of integers, (^i,-^2) 
and (^1,^2), would have been the same if, going from ^3 00 to ^3 = 0, we have traced the 
^3— dependence of the a— and /3— cycles . Within our definition of the cycles, the two pairs are 



related to each other by the SL{2, Z) transformation (|5.3|) 



(5.47) 



For qs ^ q = one finds from (|5.42|) and ( ^.34|) 



Regs 



2^3 



Re ((1 - q,/q) [In (1 - %/q) - c]) = Sni', q, - 61n2 , (5.48 



with (^1, ^2) the same as in ( ^.46| ) and (|5.47| ). It follows from this relation that Im q^/ Re ^3 ~ In q2, 
so that gs is dominated by its imaginary part at large q2- Eq. (|5.48|) defines the scaling behaviour 
of quantized gs in the vicinity of the point on the moduli space gs = —i{q2/^Y^'^- 

For g3 — >■ cxo, the leading order expression for the charge is given by (|5.45| ). One can improve 
this relation by including nonleading corrections to the periods in ( |5.42| ). In this way, we obtain 



where 



1/3 £3(2/3) 
In 



Q(n) = -(£i + 



3^3 



92 



3^3 



92 



27r |Q(n) 



2tx \Q{n) 



+ 0(g: 



E 

k=l 



Ilk e 



j7r(2fc-l)/3 



(5.49) 



(5.50) 



One verifies that the first two terms in the r.h.s. of ( |5.49| ) coincide with ( ^.13|) for = 3 and g2 
real. Eq. (|5.49| ) defines the scaling behaviour of quantized g3 in the region q^ ^ gg^^. 

By the construction, Eqs. ( |5.46| ), ( |5.48| ) and ( |5.49| ) satisfy the Whitham equation ( [5.30[ ) in 
the different regions of the (g2, gs)— space. They define the (£1, £2)— trajectories, which start at 
g2 = 1/4 and go to larger g2 = 1/4 + z/^. Example of such trajectories is shown in Figure |^. 
In the regions g2 < 1 and g2 ^ 1 the charge ga satisfies Eqs. ( p.49| ) and ( p.46| ), respectively. 
Solving the Whitham equation ( ^.301) and using (|5.49| ) as an initial condition at g2 = 1/4, one 
can reconstruct the flow of g3 in the intermediate region g2 ~ 1 along the trajectories labelled 
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Figure 8: The functions 
from ( |5.46|) using the exact eigenvalues gs. Pairs of integers 
different trajectories for the charge gs = g3(£i,£2), Eq- (|5.17| ). 



(3,0) 
(2,0) 

(1,-1) 



1/2 
12 



3.5 



—£'2 (right panel) are calculated 
'2) attached to the curves specify 



by integers (£1,^2)- Some of these trajectories go in the complex (g2, gs)— space in the vicinity of 
the point ^3 = — i(g2/3)^/^. In particular, this is the case for the (1, —1)— trajectory. The flow of 
ga around this point, shown by cross in Figure |^ (see left panel), is described by Eq. ( |5.48D . 



To verify the WKB quantization conditions, one substitutes the exact values of ga into the 
l.h.s. of ( p.46| ) and calculates the corresponding values of £1,2 for different g2 using ( ^.471) . In this 
way, for each trajectory gs = q^{v\ £1, £2) we obtain two functions that we denote as ^'i2{(l2)- Few 
examples of such functions are shown in Figure From ( |5.46|) one would expect that for large g2 



these functions should approach the same integer values £1,2 as those specifying the trajectories. 
Indeed, one finds from Figure ^ that this happens for all trajectories except the (2, 0)— trajectory. 
In the latter case, the function (i'f^\q2) approaches the value 1 instead of expected (.1 = 2. The 
reason for this is that the charge % takes anomalously small values along the (2, 0)— trajectory 
(see Figure |^ on the right). As a consequence, the nonleading WKB correction to the "action" 
function in ( p.2|) becomes important for this particular trajectory. It provides a contribution to 
the a(g3) and anils) comparable with the leading order expression ( |5.46|) and increases the value 
of i[^'\q2) improving an agreement with the exact result. 
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Figure 9: Ratio of the exact, (tz^ ■, and the asymptotic, ({^^^ ■, expressions for the charge g3 at 
small g2. Different curves correspond to the same trajectories as in Figures |^ and |[ 
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As can be seen from Figure Eq. (|5.46|) is not satisfied at small q2- To describe the Whitham 
flow of the charge gs in this region, one has to apply (|5.49|) . Comparison of the exact gs with 
the asymptotic expression ( p.49| ) for g2 < 1 is shown in Figure ^. One observes that, aside from 
the (2, 0)— trajectory, the agreement is rather good. The accuracy can be further impoved by 
including nonleading WKB corrections in ( |3.2| ). 

Thus, the WKB quantization conditions, Eqs. (|5.46| ) and ( ^.49| ), successfully describe the 
exact spectrum of the charge shown in Figure |^ for 1/4 < g2 < 1 and g2 > 1, respectively. 



6. Conclusions 

In this paper, we have developed a quasiclassical approach to solving the spectral problem for 
the noncompact SL{2, C) Heisenberg spin magnet. It allowed us to understand hidden symmetry 
properties of the energy spectrum. We also demonstrated that the energy spectrum obtained 
within this approach is in a good agreement with the exact results. 

The model represents a generalization of the well-known spin— 1/2 XXX chain to infinite- 
dimensional representation of the SL{2, C) group. Using realization of spin operators as differen- 
tial operators acting on the plane, one can map the noncompact spin magnet of length N into a 
two-dimensional completely integrable quantum-mechanical model of particles with a nearest- 
neighbour interaction. This model has appeared in high-energy QCD as describing multi-gluonic 
compound states in the multi-colour limit. Due to its complete integrability, the energy spectrum 
of the particle system is uniquely specified by the total set of quantum numbers q2, . . . ,qN- 
The latter are defined as eigenvalues of the mutually commuting integrals of motions and their 
possible values are constrained by the quantization conditions. 

Applying the methods of nonlinear WKB analysis [^, we constructed the wave function of 
the A^— particle system in the representation of the separated coordinates. To the leading order 
of the WKB expansion, the wave function is determined by the "action" function, which satisfies 
the Hamiltonian-Jacobi equations in the underlying classical model. In the classical case, the 
noncompact magnet describes the system of interacting particles moving on the two-dimensional 
plane. Solving the classical equations of motion, one finds that their collective motion describes a 
propagation of the soliton wave in the closed chain of particles with periodic boundary conditions. 
The same motion in the separated coordinates corresponds to wrapping of classical trajectories 
around the Riemann surface defined by spectral curve of the model. The charges q2, . . . ,qN 
take arbitrary complex values in the classical model and define the moduli of Ftv- 

The quantization conditions for the charges q2, ■ ■ ■ ,qN follow from the requirement for the 
wave function of the A^— particle system to be a single-valued function of the separated coor- 
dinates. To the leading order of the WKB expansion, these conditions have the form of the 
Bohr-Sommerfeld relations imposed on periodic orbits of the classical motion on the spectral 
curve Ftv- Solving the WKB quantization conditions, we demonstrated that, for fixed total 
SL{2, C) spin of the system, the eigenvalues of the integrals of motion form a lattice structure 
on the moduli space of the model. At A^ = 3 and A^ = 4 the lattices are built from equilateral 
triangles and squares, respectively. The dependence of the charges on the total SL{2, C) spin is 
governed by the Whitham equations, which were solved at A^ = 3 by making use of the modular 
properties of the elliptic curve F3. 

A novel feature of the obtained quantization conditions is that they involve both the a— and 
/3— periods on F^r. Notice that in conventional one- dimensional lattice integrable models, like the 
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SL{2,'R) Heisenberg spin magnet and the Toda chain model, the WKB quantization conditions 
involve only the a— cycles, since the /5— cycles correspond to classically forbidden zones. This 
implies that the quantization conditions for the SL{2, C) magnet are invariant under modular 
transformations of the spectral curve. As a consequence, the energy spectrum of the model 
possesses a hidden symmetry which is analogous to the 5— duality in the Yang-Mills theory ||29|| . 

In conclusion, we should mention that our consideration was restricted to the leading order 
of the WKB expansion. The obtained expressions for the energy spectrum can systematically 
improved by including nonleading WKB corrections. 



Acknowledgements 

We are most grateful to A. Gorsky and J. Kotahski for collaboration at the early stages of 
this project. We would like to thank R. Janik, I. Kogan, F. Smirnov and A. Turbiner for useful 
discussions. This work was supported in part by the grant 00-01-005-00 of the Russian Foundation 
for Fundamental Research (A.M. and S.D.), by the Sofya Kovalevskaya programme of Alexander 
von Humboldt Foundation (A.M.) and by the NATO Fellowship (A.M.). 



A Appendix: Matching conditions 



Let us require that the obtained asymptotic expressions for the Baxter blocks, Eqs. ( plj) , (|4l3|) 
and ( [4.15|) , have to verify the Wronskian condition ([4.2|) . Taking into account (|4.15|) , we obtain 
after some algebra the following expression for ( [4.2|) p| 



TT 



(A.l) 



sin(7r(s — iu)) ■ (p+{u) (</)_(«*))* — sin(7r(s + iu)) ■ (p-{u)[(p_^_{u*)y 

Together with ( [4.19| ) and ([4.21|) , this leads to the following relation for the a— functions 

sin(7r(s — iu)) ■ a+(u)(a_(n*))* — sin(7r(s + iu)) ■ a_(u)(a+(n*))* = const , (A. 2) 

which should hold for arbitrary complex u. 

Another constraint on the a— functions comes from Eqs. (|4.5| ) and (^4.1|) . It is easy to see that 
at (m = i{l — s) + e,U = —is + e) and {u = —is + e,U = i{l — s) + e) the blocks ( [4.11[ ), ( |4.13| ) 
and ( [4.15|) have poles in e generated by the first term involving functions, while the second 
term proportional to functions is suppressed as e^. As a consequence, substituting ( |4.1| ) 



into (|4.5| ) we find that (p+{u)(p+{u) — e^*° {ip+{u*)(p-^-{u*))* has to scale as ~ e as e — around 
the above two points. Applying Eqs. ( [4.19| ) and ( |4.21| ) and taking into account ( [4.22] ) we get 

a+(±2s + e)a+{±is + e 



+ 0{e''). 



(A.3) 



(a+(±2s + e)a^{±is + e))' 
As we will see in a moment (see Eq. ( [A.ll| )), this relation is exact and it holds for arbitrary real 



e. This implies in particular that S is real in (j^A_ 



So far, we have obtained two different asymptotic expressions for the function Q{u, u). One of 
them follows from the WKB expression for the wave function, Q{x/T],x/r]), Eq. ( p.l5|) . Another 

^In arriving at this relation we neglected terms suppressed by a small parameter ( [4.10|) . 
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one follows from ( [4.1|) after one replaces the Q— blocks by their expressions, Eqs. ( [4.11|) , (|4.13|) 
and (|4.15|) . We remind that the two expressions were obtained for large values of the charges g„, 
but in the different regions of parameters, x/r] 3> 1 and M=fixed, respectively. Choosing u = x/rj 
and u = X* /rj, we notice that Eqs. ( p. 15 ) and ([4.1| ) have to coincide for a; <^ 1 and m ^ 1. 

To perform the matching, we examine the behaviour of the holomorphic wave functions 
Q±{x/ri) in (|3.15|) ioi u = x/rj = fixed as — > 0. Using Eqs. (|3.12|) , and (|3.1CI| ) and choosing 
xo = 0, one finds after some algebra 



Q±{u)^ '/'exp(±zT(n)) 
where the notation was introduced for the function 

tN{u) _ 

fc=l 



N 



N 



T{u) =u In ln(AA: - m) + ^ In A/, , 



(A.4) 



(A.5) 



k=l 



with Afc being the roots of the polynomial tAr(u) defined in ( [4.16|) and ( [4.17|) . Notice that T(0) = 0. 
To obtain Q±{x/ri) one has to replace in ( [A.4| ), t]\r{u) — * i]\r{u*) = {t]\r{u))*, s —>■ s = 1 — s*, 
Afc ^ Afc = A^. Then, substituting ( |A.4| ) into (|3.15|) and making use of ( |3.21| ), one gets 



Q{u,u*) = const X |t^(M)|-^exp {-iN Arg {u^'-^)) cos (2 Re Tl 



u 



6) 



(A.6) 



,2s-l^ 



-i ln{u 



2s-l^*2s-l 



)/2. This expression should match into ( [4.1D at large u and 



where Arg (u 
u = u* . 

Let us find the large— m asymptotics of the Q— blocks, Eqs. ([4.11|) , ( [4.13|) and ( [4.15|) . To stay 
away from the poles of these blocks on the complex m— plane, we choose Re(l — s + iu) > 0. One 
gets from ( [4.20| ) (up to an inessential overall constant) 



^+(m) r^(l-s + m) 



where T(m) is given by ( [A.5| ) and 'd{u) is defined as 

d^u) = muN_ + Afcln(zAfc)+ Afcln(— iAfc), 

ImAfc<0 ImAfc>0 



(A.7) 



with N_ introduced in ( [4.23| ). We observe a striking similarity between ( |A.7|) and ([A.4|) . Identi- 
fying the r.h.s. of the first and the second relations in ( |A.7| ) as Q^{u) exp(±z'i9(u)), respectively, 
we get from ( 4.11|) , (^4.13|) and ( [4.15 ) the following expressions for the holomorphic blocks 



Q'f\u) ~ a_(iz)g+(M)e-'^(")4 
Qt\u) ~ (a_(M*))*g+(M)e-^''(")4 



sm TT s 



iu)) 



TT 



a+(M)g_(M)e'''(") 



sin(7r(s + iu)) 



(a+(M*))*g_(M)e^''("), (A.9) 



and similar expressions for the antiholomorphic blocks 



[U] ~ 



a4n)g4n)e-(-^(^*»'+ ^^^^''^'' + '''^^ a+(n)a(n)e--W^'))* 

TT 



Q'f'\u) ~ (a_(M 



*))*g_(^)e-W"*»' '^^^ (a+(r))*g, (^)e-w^*))' 

TT 
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(A.IO) 



The functions Qo^{u) and Q±{u) approximate the exact solutions to the holomorphic Baxter 
equation ( |2.11| ) in the different regions, m ^ 1 and u ~ 1, respectively. Eq. ( [A .91 ) sews these two 
sets of functions in the intermediate region of u. 

Let us substitute Eqs. ( |A.9|) and (|A.10|) into (4J.) and compare the resulting expression for 
Q{u,u) with ( |3.15| ) at m = x/r] and u = u*. One finds that Q{u,u*) involves four different 
combinations of the Q— functions, Q±Q-i- and Q±Q::^, whereas the r.h.s. of (|3.15|) contains only 
the diagonal terms. To cancel the off-diagonal terms ~ Q±Q^, one has to require that 



(a+(n)a+(n*))* (a_(n)a_(M*))* 



The coefficients c± in front of the diagonal terms Q±Q^ are given by 
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— e 

IT 



2iRei9{u)-i(5 



sin(7r(s + iu*))a^{u)a+{u*) 



sm{7T {s -iu*)){a+{u)a-{u*))* , (A.12) 



(A.ll) 



and c_ = —c*^- Eq. ( |A.12| ) can be further simplified. Replacing 5 by its expressions ( |A.11| ) and 
making use of the Wronskian relation ([A.2|) , one gets (up to an overall normalization factor) 



-2iRe-d{u)-i5 O'-ju) 



(A.13) 



Here, in the last relation we applied the identity [sin(7r(s — iu))]* = sin(7r(l — s + iu*)) = 
(— 1)"^"'"'^" sin(7r(s — iu)) with n„ = i[u — u*) integer in virtue of ( |2.10| ). Finally, the resulting 
expression for Q{u, u*) matches (|A.6|) provided that the coefficients c±, Eq (|A.12|) , do not depend 
on u. In addition, these coefficients satisfy the relation ( |3.21| ), c+/c- = —c+/c\ = exp (—2^9), 
which leads together with (|A.13| ) to 



=Fie 



a_(u] 



{a4u*)r 



I e 



2iKe^{u)+i5-i@ 



(A.14) 



We recall that these relations hold for u taking the same values as separated coordinates (|2.10|) , 
that is u = Vu — with Uu integer and real. For a given set of the integrals of motion, 

the phases 'Re'diu) and 9 entering (|A.14|) are uniquely by Eqs. ( |A.8| ) and ( p.24| ). In contrast, 
the phase 5 depends on the normalization of the blocks. Substituting a± (u) e*''/^ a± (u) and 



a±{u) e*"^/^ a±{u) in ([A.14|) , one can put 5 = in Eqs. pTTp . 



Eqs. (|A.14|) and (|A.11|) fix the phases of the a— functions but not their absolute values. 
Nevertheless, as shown in Appendix |B|, this data becomes enough to construct the eigenvalues of 
the Baxter operator Q{u,u*) and to calculate the energy spectrum of the model. 



B Appendix: Calculation of the energy spectrum 

In this Appendix we obtain the WKB expressions for the energy and quasimomentum, Eqs. ( 4.25|) 
and ( |4.24| ), respectively. 

To begin with, we substitute u = ±is + e and u = ±is + e into asymptotic expressions for the 
Q— blocks, Eqs. ( |4.11| ) and ( [4.13| ), and examine the limit e — 0. One finds that the F— functions 
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in the denominator suppress the contribution of one of the terms in the r.h.s. of ( [4 .111) and ([4.13|) . 
This leads to the following expression for the quasimomentum 



a+hs)a^ns. 
9n = iln ._' — + I In 



(p+{is){^-{is-i))* _n{a),n(b) 



a+{—is)a^{—is) (p-{—is){Lp^{i — is)) 



'N ~ "N y 



(B.l) 



where, for convenience, we split the expression into two pieces. One applies (^4.201) and uses the 
asymptotic behaviour of the F— functions at large arguments to get 



e 



(b) 
N 



-4 Re 



^ Afcln(zAfc)+ ^ Xk\n{-i\k) 

,ImAfc<0 ImAfc>0 



(B.2) 



The calculation of 6^^^ goes as follows. One substitutes u = ±is into ( |A.2[ ) and uses anti- 
periodicity of the function a+{u), Eq. (|4.22|) , to obtain a+{is)(a^{is))* = a_(— 2s)(a+(— is))*. 



This relation allows us to rewrite 6^^^ as 



(a_(— is))* a+{is) 



Uj^ — bin- I ■ w'' 



-26 + 1 In 



(a_(— is))* (a+(is))* 
a_(— is) a+(is) 



(B.3) 



where in the second relation we used Eq. ( |A.11| ). Then, one takes into account (|A.14| ) and (|A.8|) 

to get 



9 



(a) 
N 



-20 + 4 Re 



^ Afcln(iAfe)+ ^ Afeln(-iAfc) 

.ImAfe<0 ImAfc>0 



(mod 27r) . (B.4) 



Finally, combining together ( |B.2|) and ( |B.4| ) we obtain ( [4.24| ). 

The calculation of the energy ( [4.6|) goes along the same lines. One substitutes ( [1.1 1| ) and 
( [4.13| ) into ( |4.6| ) and expresses the result in the following form 



E 



N 



-2 Im In 



+ A^iv ^ - 2 Im In 



a+(zsj 



[^+(z(l - s))]v J V a+[-is 

Here, the notation was introduced for an additive constant 

AEm = -2N Re[ip{l - 2s) + V^(l - 2s)] +6^ = -4Ar^(l) 



7T.(a) I 7T.(f') 



(B.5) 



(B.6) 



with En defined in ( [4.6] ) and Re'?/'(1 — 2s) = Re-?/'(2s — 1) in virtue of s = 1 — s*. To evaluate 
E^^ we apply ( [A.3| ) to get 



E 



(b) 
N 



Im ( In 



a+(is) (a+(— is))* 
(a+(is))* a+(— is) 



Im ( In 



a+(is)) (a+(— is)) 
a+(is) a+(— is) 



0, (B.7) 



where the last relation follows from (|A.14|) and (|A.8| ). Finally, one rewrites E^^^ as 

= -2 Im ( In [^+(zs) ^+(z(l - s))] )' - 4iV^(l) , 
takes into account ( [4.20| ) and arrives at ( |4.25| ). 



(B.8) 
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According to (|4.24|) , the quasimomentum is given by 



On - 



— Re 

TT 



dx Sq{x) 



(B.9) 



where the "action" differential was defined in ( p.GD H 

NtN^x) — xt'j^{x) 



dx Sq{x) 



y{x) 



dx . 



(B.IO) 



In Eq. ( [B.9|) the integration contour, 7p-p+, goes on the Riemann surface (|1.10|) from the point 
located above x = on the lower sheet to the point Pg^ above x = on the upper sheet and 

does not intersect the cycles and {k = 1, N — 2) as shown in Figure |l|. 
Let us consider the following integral 



dx S'q{x) 



7(o"fc) 



^^ Ntj^{x) - xt'^{x) 



(B.ll) 



Here the integration contour 7(0"^) goes from the point Pf^ above x = on the upper sheet of 
Fat to the branching point ak, encircles it and goes to the point P^ on the lower sheet of F^v. 
Notice that 7(0"^) is different from 7p-p+, but the two contours are related to each other through 
( p:22|) . It follows from (^) that i = N/TiReh (mod N) leading to 

In general, Ik is a complicated function of the integrals of motion q. The integral in ( [B.ll|) 
can be easily evaluated for q satisfying ( |5.7| ). In that case, the spectral curve takes the form 
( |5.8|) and the branching points are given by (|5.1CI| ). Replacing the integration variable in ( [B.ll|) 
as X ^ unx and taking into account the hierarchy (|5.9|) one gets for j = 1, . . . ,N 



Ij = 2u 



N 



.i7r(2j-l)/JV 



Ndx 



1 



X 



N 



-^xd^ )p^_2(x) + C(e2^ 



(B.12) 



7y 9 

U2X + ... + Un~iX 



where Pn-2{x) 
Straightforward calculation leads to 



1 + x^ 

C(e) and we neglected terms quadratic in Pn-2{x). 



2u 



N 



J-K{2j-l)/N - 



N-l 



n=2 



(B.13) 



with B(x,?/) = F(x)F(y)/F(x + y) being the Euler function. This expression has the following 
properties 



Ij+N 



N 



3/2 N 



(B.14) 



The r.h.s. of ( |B.13| ) takes the form of a discrete Fourier transformation from the coordinate (j) 
to the momentum {n) representation. This allows one to get the moduli and m„ as inverse 
Fourier transformation of Ij. Then, taking into account that Re/^ = Ti{nk + i/N), Eq. ([5. 51), one 



arrives at (|5.11 



''Here we neglected the last term in the r.h.s. of (3^) since it does not contribute to (B.9). 
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C Appendix: Calculation of the quasimomentum 

Let us demonstrate that for the integrals of motions q satisfying the quantization conditions 
( p.23| ), the parameter £ takes strictly integer values in ( [B.9| ). To begin with, we define on the 
Riemann surface ( |1.10|) the set of normalized differentials of the first kind, Uk, and the third kind 
(dipole differentials), Qoc, 

Ukj , ^oo = '2qr + y.Uj , C.l 

3=1 y y ~t y 

with the expansion coefficients Ukj and Uj fixed by the normalization conditions 

tuk = 2715 jk , (p Vl^ = Q. (C.2) 



The differential f2oo has a pair of poles located at the points (-Poo -^oo) above x = oo on the upper 
and lower sheets of Ftv and the residue at these poles resp± VL^o = ±1- The action differential 
( P-lOp can be decomposed over the set of the differentials (|C.1D as 



N-2 



dx S'o(x) = ql^^ -n^ + ^ak-ujk- (C.3) 



k=l 



Here, the coefficient in front of floo is fixed by the asymptotic behaviour of dxS'^lx), Eq. ( p.lO|) , 
at infinity Sq{x — > P^) ~ ±ql^'^/x. In the second term, the notation was introduced for the 
a— periods of the action differential 

°'k = ^ i dxS'oix) , Reafc = ^£2k-i , (C.4) 

whose values are fixed by the quantization conditions (|3.23|) . 

Substituting ( |(J.3| ) into ( |B.9| ) and applying the well-known identities between the contour 
integrals on the Riemann surface p2| 



p+ /. p+ « p+ 



^oo= ^0, / ujk = -i f VLq , (C.5) 

'P^ JPq Jh 



one can express £ in terms of the integrals of the dipole differential VLq 



TT 



„p+ N-2 „ 

q^^ / ^]o - ^ Ofc i fio 

JP^ k=i JPk 



(C.6) 



The differential VLq has a pair of poles located at the points on Fat above x = and is 
normalized as £j, f^o = 0. This differential plays a special role in our analysis as it can be 
expressed in terms of the quasimomentum 

= -^dxp'ix) , e^(^) +e"^(^') = ^ . (C.7) 
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Notice that p{x) is a multi-valued function on the Riemann surface ( |1.10|) such that 
with (Tfc being the edges of the cuts, t%{ak) = ^crl^ ■ As a consequence, 

r.P+ 



Tfl 

Qq = 2TTi — 



with m and integer. Finally, we find from ( |C.6| ) 



fin = 27Ti—^ 



N 



2iVRe 



.m 1/2 



Ar-2 



1/2 , \^ 

'i2 +2^ 



ruk 



k=l 



N-2 



N-2 



nm 



+ 2 nik Re 



Qk = nm 



k=l 



k=l 



rrikhk-i 



(C.9) 



(C.IO) 



Here, in the second relation we took into account that q^^"^ = i{h — 1/2) + 0{{h — l/2)~^) = 
in/2—p , with h = (l+n)/2+iz/. Notice that the quasimomentum depends only on the a— periods 
of the action differential. 
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